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Abstract 



The standard framework within which cosmological measurements are con- 
fronted and interpreted nowadays, called A Cold Dark Matter, presents a 
Universe dominated by unknown forms of energy and matter. 

My Thesis is devoted to investigate the distribution of dark matter in 
galaxies and addresses the fact that the local universe-the small objects that 
orbit galaxies and the galaxy cores-turns out to be a marvelous laboratory for 
examining the nature of dark matter and the fundamental physics involved 
in the structure formation and evolution. 

I develop tests, based on mass modeling of rotation curves, for the valida- 
tion of dark matter models on galactic scales. These tests have been applied 
in analyzing the phenomenology of the cusp vs core controversy, and the 
phenomenon of non-Keplerian rotation curves as modification of the laws of 
gravity. I further investigate the properties and scaling laws of dark matter 
halos. 

My conclusion is that galactic observations provide strong imprints on 
the nature of dark matter. 
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Preface 



Cosmological observations provide compelling evidence that about 95% of 
the content of the Universe resides in two unknown forms of energy that we 
call dark matter (DM) and dark energy: the first residing in bound objects 
as non-luminous matter, the latter in the form of a zero-point energy that 
pervades the whole Universe [I]. The DM is thought to be composed of cold, 
neutral, weakly interacting particles, beyond those existing in the Standard 
Model of Particle Physics, and not yet detected in accelerators or in ded- 
icated direct and indirect searches. In the standard A Cold Dark Matter 
(ACDM) scenario primordial density fluctuations are generated during an 
inflationary period and become the seeds of the bottom-up structure forma- 
tion model. This scenario successfully describes the accelerated expansion 
of the Universe, the observed temperature fluctuations in the cosmic mi- 
crowave background radiation, the large scale matter distribution, and the 
main aspects of the formation and the evolution of virialized cosmological 
objects p]. 

Despite these important achievements, at galactic scales of about 10 kpc, 
where today most of the mass is located, the ACDM model meets with se- 
vere difficulties in explaining the observed distribution of the invisible matter 
around the luminous one. In fact, on the one hand, N-body simulations per- 
formed in this scenario, unambiguously predict that every halo in the Uni- 
verse hosting and surrounding a galaxy, must have a very specific density 
profile. This features a well pronounced central cusp, obeying to the well 
known Navarro, Prenk & White (NFW) profile [3j. On the observational side 
instead, high-resolution rotation curves (RCs) show that the actual distri- 
bution of DM is much shallower than the above, and it presents a density 
profile with a nearly constant density core [H El [H [3, [TOl [H] that is well 
represented by a Burkert profile 

The cusp vs core controversy, together with other present main failures 
of the the standard collisionless particle paradigm, such as the issue on the 
number of sub-halos [12], has far-reaching consequences in the researches on 
the nature of DM and unveil the elusive knowledge on it. 

My aim is to show how a systematic comparison of cosmological and 
particle physics models with galactic observations provides strong bounds 
on the properties of DM. The outline of this Thesis is then as follows. In 
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the Introduction, the basis of the ACDM scenario is summarized, presenting 
fundamental physical ingredients and its important predictions. In Chapter 
2 a picture of the fundamental properties of DM as non-relativistic particles 
is given, as well as their devoted searches. There are in fact hints that the 
phenomenon of the missing mass is linked to a new high energy phenomenol- 
ogy not included in the standard model of particle physics and foreseeing the 
existence of new elementary particles with a mass roughly above the hundred 
GeV scale. 

Chapter 3 is devoted to a discussion on the mass modeling of RCs as 
major tools for investigating the distribution and nature of DM in galaxies. 
A deep understanding of the mass models is a fundamental ingredient for 
a future discovery of the DM particles or alternatives theories, in what it 
provides the link between the microphysics phenomenology and the obser- 
vations. Moreover this chapter gives an exhaustive discussion on the current 
status of the cusp vs core controversy with extensive references to relevant 
literature. 

In Chapter 4 basic ideas behind the most popular alternatives to DM, 
MOND and f(R) theories of gravity, are summarized. In Chapter 5 I develop 
a test for analyzing the models which aims at solving the cusp vs core con- 
troversy by resorting to the best available galaxy kinematics. I apply this 
test to an example for such models. In Chapter 6 instead I develop a test 
for analyzing the models which aims at solving the phenomenon of the RCs 
by resorting to modifications of the laws of gravity. I focus on f(R) theories 
of gravity. This work represents a step forward on the issue in what for the 
first time a complete analysis with a devised RC sample has been performed. 

I have also further investigated the properties of DM halos. Kinematic 
observations of the dwarf spheroidal (dSphs) satellites of the Milky Way are 
revealing hints about the structure of DM halos. I investigate whether the 
extrapolation of the scaling relations of brighter galaxies to the low end of the 
galaxy luminosity regime is consistent with the observed internal kinematics 
of dSphs. In Chapters 7 and 8 I discuss the implication of such relations for 
the comprehension of the nature of DM. Finally I conclude in Chapter 9. 
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Chapter 1 

Introduction 



More than eighty years ago E. Hubble estabhshed the expansion of the Uni- 
verse with his pioneering observations of galaxies. Since then galaxies have 
been fundamental tools for understanding the structure and evolution of our 
Universe. Today they are crucial laboratories where microphysics phenom- 
ena, up to now not detected by particle physics experiments, emerge with 
unprecedented clarity. In particular the great improvement in quality and 
quantity of the measurements of galaxy kinematics spanning a large range 
in luminosity, has provided precise tests for evaluating theories both of cos- 
mological and particle physics relevance. 

The study of the micro and macro cosmo today produced the ACDM 
scenario (not yet a theory!) which allows the study of the formation and 
evolution of cosmic structure from first principles, and embraces cosmological 
theories (Big Bang and Inflation), particle physics models (the standard 
model and extensions) and astrophysical models and observations. 

The fact that we need a mysterious new form of matter having a domi- 
nant role in structure formation and evolution represents for the first time a 
demonstration from the cosmological side that the standard model of particle 
physics needs a deep extension. 

In this Chapter I first give a brief introduction to modern cosmology 
with reference to the latest precision measurements of its most important 
parameters. The first paragraph provides the cosmological basis for the 
ACDM paradigm. Then I introduce the theory of structure formation and 
the growth of perturbations in the primordial Universe. I finally describe 
the particle physics basis aspects relevant for this paradigm. 

1.1 Observations of the Large Scale Structure of 
the Universe 

Within the current ACDM paradigm of structure formation and evolution, 
cosmology provides the initial and boundary conditions that together with 
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astrophysical models allow to make definite predictions about the visible 
Universe. The systematic comparison between these predictions and the 
astrophysical observations are fundamental tests of any cosmological model. 

Our modern theory of the universe, started with the work of Einstein 
and Friedman in the 1920s, is based on the Einstein's theory of space-time 
developed few years before. It starts from the assumption of homogeneous 
and isotropic universe at large scales, described by the Friedmann-Lemaitre- 
Robertson- Walker (FLRW) metric: 

ds'^ = dt^ - R^(t) \ — r + de^ + sir? 6 #H , (1.1) 

[ 1 — A; J 

where {t, r, 6, (f>) are co-moving coordinates, i2(t) is the cosmic scale factor, 
and fc is a curvature parameter which can be chosen to be +1,0 or —1 for 
positive, flat or negative curvature respectively. For a test particle moving 
freely in such a metric the geodesic equation reduces to: 

\p\ dt R dt' ^ ' 

where p is the particle momentum. This equation shows that the relativistic 
momentum is red-shifted by an amount z = — 1 as the scale factor 

expands. 

In the Big Bang model the scale factor evolves over time and its evolution 
is related to the energy density by the two Friedmann's equations: 

R k SttG 

« + « = 

and 

R AttG , , , , 

= ^(p + 3p), (1.4) 

where G is the gravitational constant and p and p are the energy density and 
pressure of the universe. In a Newtonian interpretation the first equation is 
the energy balance in a central force problem, while the second one is the 
analogous of the Newton law F = ma. It is possible to define the critical 
density of the universe as: 

= 8^' ^'-^^ 
which corresponds to the density of a flat Universe. 

The energy density and pressure in general receive contributions from 
several kind of sources like photons, baryons, DM and several others: 

p = + Pb + PDM + ■■■ P = P-r +Pb+ PDM + ■■■ ■ (1-6) 



From general thermodynamic reasoning however all these sources respect 
a general relation between density and pressure: p = wp, where w = 1/3 
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is valid for an ultra-relativistic fluid (radiation), tt; ~ is valid for non 
relativistic species (matter) and w = —1 is valid for vacuum energy. This 
implies that the energy density of radiation scales with the expansion of the 
universe as p-y oc while for the non relativistic components and for 

vacuum energy we have respectively: pcDM and px oc R^. It follows 

that the early Universe was dominated by the radiation energy density while 
at later stage it became matter dominated, with estimates of the redshift 
of the transition epoch of z ~ 10^ when it was about i ~ 5 • 10^ years 
old. The time of radiation-matter equality is of fundamental importance for 
the understanding of the formations of the structure of the Universe as it 
represents the moment when the primordial density fluctuations start to have 
a significant growth. In the last stages vacuum energy however dominates. 

From the 1970s the FLRW cosmology is rooted in three observational ev- 
idences: the expansion of the Universe, discovered in the 1930s by E.Hubble 
observing the recession of galaxies as a function of their distance. The sec- 
ond evidence is the Primordial Nucleosynthesis, pioneered in the 1940s by 
G.Gamow. Finally the Cosmic Microwave Background (CMB), which is 
the fossil radiation of the primordial universe discovered in the 1960s by 
A.A.Penzias and R.W.Wilson and today has a temperature of Tq ~ 2.7 K. 

The Hubble law is described by the equation V = Hod, where V is the 
recession velocity, d is the galaxy distance and Hq ~ 71km s^^ Mpc^^ is 
the Hubble constant. The Hubble constant is linked to the scale factor R by 
a Taylor expansion: 

= 1 + (t - to) - IqoHi {t - tof + ... (1.7) 

where qo is the so-called deceleration parameter linked to the second deriva- 
tive of the scale factor. 

The Primordial Nucleosynthesis explains the relative abundances of light 
elements (Hydrogen, Deuterium, Helium-3, Helium-4 and Lithium-7, see 
Fig. Il.lj) produced during the first 20 minutes of the Universe. The pre- 
diction depends on one free parameter: the baryon-to-photon ration r/ = 
273-10~^^Qbh^ where is the ratio of the baryon density to the critical 
density and h defined such that Hq = 100/ikm s~^ Mpc~^. Measurements 
give 4.7 • 10~^° < Tj < 6.5 • 10"^" ^3], giving a precise measurement of the 
baryonic content of the Universe. Moreover these measurements yield an 
Helium-4 mass fraction Yp ~ 0.25. This mass fraction is of great importance 
in developing the mass models of gaseous disks in spirals. 

Despite the success of the FLRW cosmological model, the Universe is 
clearly neither homogeneous nor uniform on scales smaller than about 100 Mpc 
The modern trend in cosmology, both theoretical and observational, is to try 
to understand the formation and evolution of the inhomogeneities of cosmo- 
logical relevance. What follows is a brief description of the most important 
observations on the very large scales. 
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Figure 1.1: Solid lines: predictions of light elements abundances from BBN. 
Shaded areas: best fit for the fundamental baryon-to-photon parameter. 
Rectangles with dashed contours: measurements of light elements abun- 
dances [Hj. 



After the discovery by G. Smooth and J.Mather with the COBE mission 
of fluctuations in the CMB temperature of the order 6T/T ~ 10^^ (after the 
subtraction of the variation due to earth motion of order 6T/T ~ 10~^), these 
measurements played a major role in the development of the cosmological 
model and of the ACDM paradigm of structure formation. The importance 
relies on the fact that the fluctuations in the CMB temperature reflects the 
fluctuations of the matter density at the time when the primordial plasma 
became neutral, at a temperature T ~ 3000 K and redshift z ~ 1100: 

— ~ const X (1.8) 

where the constant depends of the kind of matter considered. Up to the 
recombination epoch the temperature fluctuations evolve under the influence 
of sound waves propagating in the hot plasma. 

Fig. 11.21 shows the more recent CMB measurements from WMAP. From 
this figure fluctuations with an angular size of about half a degree are clearly 
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Figure 1.2: The full sky 5-years WMAP image of the CMB temperature 
fluctuations after dipole subtraction and foreground reduction |15| . 

6000 I 1 1 1 1 1 — I 




10 100 500 1000 1500 
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Figure 1.3: The WMAP 5-year power spectrum along with recent results 
from the ACBAR, Boomerang, and CBI experiments [15]. The red curve is 
the best-fit ACDM model to the WMAP data. 

visible. A quantitative analysis performed on the multipole decomposition 
of the 2-point correlation function (see Fig. II. 3|) shows that multipoles with 
I < 10^ corresponds sound waves with periods bigger than the age of the 
Universe at decoupling. The multipoles with 10^ < / < 10^ show clearly 
the oscillations of the sound waves with period short enough to undergo at 
least one oscillation before the decoupling. The position of the first peak in 
this region is sensitive to the flatness of the Universe while the ratio of the 
height of the even peaks with respect to the odd ones gives a measurement 
of the ratio between DM and Baryon content of the Universe. Multipoles 
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Figure 1.4: Two dimensional distributions of galaxies within 1.25° of the 
Equatorial plane [16| . 



with / > 10 are suppressed due to the fact that the recombination did 
not happened instantaneously, but the last scattering surface had a finite 
thickness. 

The second pillar of the modern observational cosmology relies on the 
measurements of the galaxy distribution over large portions of the visible 
Universe (see Fig. 11.41) . One of the main challenges of any theory of struc- 
ture formation and evolution is to explain how the tiny fluctuations in the 
baryon density measured by the CMB evolved under the influence of gravity 
up to the stage visible today within the known age of the Universe (this is 
precisely one of the main success of the ACDM paradigm). From the Fourier 
analysis of the two-point galaxy correlation function (see Fig 1 1.51) it has been 
possible recently to measure the imprint of the primordial sound waves (in 
this context called Baryon Acoustic Oscillations, BAO) in the visible Uni- 
verse. The primordial fluctuations start to oscillate due to the interplay of 
the pressure of the hot plasma and the attraction of gravity as soon as their 
size is below the horizon of a given epoch. At the time of the baryon-photon 
decoupling the plasma becomes neutral and pressure drops arresting the os- 
cillations and leaving only gravity as dominant force. The imprints of the 
primordial oscillations however is still visible in the large scale matter dis- 
tribution: galaxies in fact are encountered more often in the large overdense 
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Figure 1.5: Comparison of current power spectrum with observations from 
CMB, galaxy surveys, cluster, lensing and Lya forest [16] . 

regions than in the depleted ones. 

The third fundamental observation of modern cosmology is measurement 
of the distance of Type la supernovae which allowed the discovery of the 
accelerated expansion of the Universe. 

Other important measurements include the matter power spectrum as 
measured from Lyman Alpha absorbers and the cosmic shear (cosmological 
weak lensing) . The Lyman Alpha forest in particular is the sum of absorp- 
tion lines arising from the neutral hydrogen Lyman Alpha transitions and is 
visible in the spectra of distant objects (see ll.6|l . These absorption lines are 
due to clouds of neutral hydrogen which the emitted light encounters while 
traveling to earth. Their amplitude and position depend on the matter den- 
sity as a function of the redshift and hence is a good probe of the matter 
power spectrum. 

The combination of the CMB, the BAO and the SN la data are well fitted 
by the ACDM cosmological model, whose free parameters together with the 
best-fit values are shown in table 11.11 The same data imposes important 
constraints on models with extended sets of free parameters. The emerging 
global picture is a universe with an energy density dominated today by the 
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Figure 1.6: Spectra of two Quasars. Top: near one. Bottom: a distant one 
featuring Lyman alpha absorption. 

vacuum energy (for about 72%). Important contributions arise also from 
CDM (about 22%) and Baryons (about 4.6%), while for Neutrinos an upper 
limit of fijy < 0.026 is obtained. 

1.2 Structure Formation 

The paradigm for structure formation finds its roots in the pioneering work 
by Peebles (see e.g. [E]), who developed the hierarchical clustering theory. 
In this scenario, structure builds up through the aggregation of nonlinear ob- 
jects into larger and larger units. In the current CDM model, the build-up 
of structures is governed by the dark dissipationless component, that evolves 
under gravity from an initially gaussian distribution of primordial pertur- 
bations; small fluctuations first, and then larger and larger ones, become 
nonlinear and collapse when self-gravity dominates their dynamics, to form 
virialized, gravitationally bound systems. As larger perturbations collapse, 
the smaller objects embedded in them cluster to form more complex pat- 
terns. In the meanwhile, the DM provides the potential wells within which 
the gas cools and forms galaxies under dissipative collapse. 

The development of the proper description of the origin of the structures 
(e.g. [18j) needs two basic inputs: the initial values of the matter density 
fluctuations and a proper model for their evolution. 

The fluctuations are described by introducing the density contrast: 



6{x) 



6p{x) _ p{x) - p 



(1.9) 



P 



P 
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parameter 


ii7"i\/rAp 1 "R A O ^ Q1\T 


short description 


1 nno, /i^ 


z.z / It u.uoy 


nU.IlQ16Q 1111160 1116 DdiyOIl Q6Ilbliy 




n 1 1 4'^ -1- n 00*^4 


V_>U1Q UdlK IVldjlLcl Qclloliy 




n 791 4- n ni ^ 

U. / Zi ± U.UiD 


Dark Energy density 




u.ybu_ooi3 


Scalar Spectral index'^ 


T 


U.Uo4 ih U.UiD 


Reionization optical depth 


A 2 //^eA 


(n /I c 7+0.092 \ ^ in-9 


Amplitude of curvature perturbations^ 


0-8 


U.817 ± U.UzD 


Galaxy nuctuation amplitude 


-no 


7n 1 -1- 1 "XlCm la l'\l\-r\n 

(U.i It i.oivni/ s/ ivipc 


Hubble constant 


^reion. 


10.4 lb 1.4 


Redshift of reionization epoch 


to 


13.73 lb 0.12 Gyr 


Age of the universe 




0.0462 lb 0.015 






0.233 lb 0.013 






0.1369 lb 0.0037 


Matter density 



a) estimated at /cq = 0.002/Mpc 



Table 1.1: Summary of the cosmological parameters of the ACDM model 
and corresponding to 68% intervals from [15]. 



where p{x) is the local matter density and p its average, or alternatively, by 
the Fourier coefficients of the density contrast defined by: 



h = y~ \ S{x)exp{ik ■ x)d'^x, (1.10) 

J Vol 

where V is an appropriate normalization volume, and k is its comoving 
wavenumber. Accordingly the physical wavenumber is kphyg = k/R{t) and 
the physical wavelength is then \phys = R{t)X = 27rR{t)/k. The density 
perturbations are also characterized by the mass within a sphere of radius 
A/ 2 given by: 



TT . 



M ^ -KhysPm ^ 1-5 X lO^Mf0(n„/i^)A^p^, (1.11) 

implying that a galactic mass perturbation corresponds to a scale of a Mpc. 

The primordial fluctuations are generated randomly according to a dis- 
tribution which is considered as a power law spectrum: 5k — AVk"', where 
A is its characteristic amplitude. It is useful to introduce the root mean 
squared density fluctuation as: 

^ = {6{x)6{x))'^^ (1-12) 
which, taking the Fourier transform reduces to: 

P J Jo 27r^ k 
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It is customary to define P{k) = \Sk\^ as the power spectrum. 

Fluctuations are normally divided in two classes: curvature (or adiabatic) 
and isocurvature (or isothermal). The former are authentic fluctuations in 
the matter density while the latter are fluctuations in the matter composi- 
tion (e.g. variation in the fraction of baryons) which results in variations in 
the local equation of state. The difference between the two types however 
is relevant only on scales larger than the horizon as on smaller scales mi- 
crophysics process can transform isothermal in adiabatic fluctuations (and 
viceversa). In the following only curvature fluctuations will be considered. 

To start the study of the linear description of the perturbation evolutions 
the simple case of fluctuations in a non expanding universe will be considered 
first. This analysis allows the introduction of a fundamental quantity called 
Jeans Length. This simple analysis presents however some inconsistencies 
which can be eliminated in a more complex and rigorous model. 

In Eulerian coordinates of a non expanding Universe the equations de- 
scribing matter and momentum conservation and the Poisson equations take 
respectively the following forms: 

^+Po^-vi = 0, (1.14) 



On 2Y3 

dt ^""^ po 



+ vi-^ +V(l) = 0, (1.15) 

^0 

VVi = 47rGpi, (1.16) 



where p is the matter density, p and v its local pressure and velocity re- 
spectively, and gravitational potential. The subscript indicates the ho- 
mogeneous case (i.e. po = const, po = const, (po = const, vq = 0) and the 
subscript 1 the small perturbations (i.e. p = po+ pi, p = Po+Pi, 4' = ^^'o+^i, 
V = vo + vi)- Us = ^g^^ — ^ sound speed (in adiabatic conditions). 

The equations of the perturbations can be combined in a second order 
differential equation of the form: 

^-^;2vV = 47rGpoPi. (1.17) 

Assuming solutions of the form pi{f,t) = Ae^~'^'^"^'^'^^^^ po, the dispersion 
relation is obtained: u>'^ = Vgk^ — AttGpq, with k = |^|. 
Defining the critical Jeans wavenumber as: 



^ 47rGpo y/^ 



%=Rr) , (1-18) 



it is clear that solutions with wavenumber less than kj are unstable (either 
exponentially growing or decaying) while solutions with bigger wavenumbers 
have oscillatory behavior. 
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Considering the unperturbed solutions for the matter density, matter 
velocity and gravitational potential, in an expanding Universe, a second 
order differential equation for the Fourier transform of the density contrast 
is obtained: 

4 + 2^4 + f ^ - 4^Gpo Vfe = 0. (1.19) 



In a flat matter-dominated model, the solution of this equation for the un- 
stable (k « kj) growing (S+^k) or decaying (S-^k) mode takes the form: 



5+,k{t) = S+,k{ti)\^-j , (^_,fc (t) = (t,) l^-J , (1-20) 

where 6+^k{ti) and are the initial values at a chosen reference time ij. 

The exponential evolution obtained in a non expanding Universe becomes a 
power law evolution in an expanding Universe. A realistic treatment of the 
evolution of the perturbations however must consider the dynamics of several 
fluids, each with a different equation of state. Moreover the full treatment 
of the general relativity formalism must be taken into account. 

In the ACDM model the evolution of the linear power spectrum is con- 
structed as: 

k^P{k, z) _ ^ ^^_9 / ' X D\k, z)T\k) ( ^ 



27r2 ■ yhHlQ.^) ' \kwMAP, 

(1.21) 

where D{k, z) and T{k) is the linear growth rate an the matter transfer 
function (e.g., [Ej). The model with spectral index ~ 0.96 fits the data, 
indication an almost-free power spectrum0. Notice that the requirement of 
hierarchical clustering, that small objects form first, is ensured if P{k,z) is 
a decreasing function of mass, or correspondingly, an increasing function of 
the spatial wavenumber k. Using the fitting functions for D{k,z) and T{k) 
as found in (SHI I plot in Fig ll.7l the current power spectrum. 

The rms amplitude of mass fiuctuations inside a top hat spherically sym- 
metric window of radius R is 

/•oo dk k^ / 3 \ ^ 

a'^iR,z) = J ——P{k,z)ij—^{smkR-kRcoskR)] , (1.22) 

where Mr = ^pqR^ is the mass enclosed in the window and at the mean 
density po of the Universe (see Fig ll.Tp . 

The mass function can then be defined as 

Mdn(M,z) 



^Inflationary models favour a running spectral index, ns{k) = dlnP{k)/dlnk. 
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Figure 1.7: Prom left to right: current linear spectrum, rms of mass fluctu- 
ations and mass function. 

where n{M, z) is the abundance of halos with mass less than M at redshift 
z. In Fig. 11.71 1 plot the mass function at z = using the formula from the 
simulations of [H], f{M) = 0.315exp(-|/nCT~^ + 0.61|3-^). 

Correcting the linear prediction for the nonlinear dynamics when the 
density contrast grows above unity (important at small scales) the agreement 
with the observations is remarkable as shown in Fig. 11.51 

The study of hot DM models have been well motivated by the presence 
of neutrinos in the primordial universe and by the difficulties presented by 
CDM. A hot DM species however is characterized by a typical length (called 
free-streaming length) which is of order of several Mpc. The hot DM compo- 
nent would dump the fluctuations on scales smaller then the free-streaming, 
a disfavoured situation. 



1.3 Summary of thermal history of the Universe 

A major achievement in the development of cosmology has been to show how 
the behavior of the Universe on the large scale is dictated in a good part by 
microphysics (see e.g. [IH])- The microphysical laws are incorporated in 
the study of cosmology and structure formation in a statistical sense by the 
use the phase space distribution function fi{p^-,x^) of the various species of 
particles (i) considered. 

The evolution of the distribution functions is determined by the Boltz- 
mann equation in its general relativistic form: 

K|^-apVi^ = e|/l. (1.24) 

where V^^ is the Christoffel symbol and C [/] represent the collision term. 

The number density, energy density and pressure of particle specie can 
be obtained by integrating the distributions (using natural units and the 
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and 



relation E'^ = p'^ + m^): 

<^') = ji^ j fi^',^)d'p, (1-25) 
Pixn = j E{-^)f{x^,^)d^p, (1.26) 

P{xn = ji^ J ^^fi^'^P)d'p, (1-27) 

where g is the number of internal degrees of freedom of the particle. 

In an homogeneous and isotropic universe the distribution function is a 
function of only energy and time f{E,t), and in the equilibrium condition 
takes the form: 

f{E,t) = ^ , (1.28) 

e T rb 1 

where the temperature T and the chemical potential /x are functions of time 
and the sign is positive for fermions and negative for bosons. 

Particles of specie tjj are kept in thermal equilibrium by interaction pro- 
cesses, which for concreteness can be considered of the form ipip <^ XX, 
where X is any kind of final state particle (the elastic scattering is a simple 
example of process of this type). In this case the Boltzmann equation in a 
FLRW metric reduces to: 

^ + SH{t)n = -{av){n'-n%), (1.29) 

where H{t) = R{t)/R{t), a is the cross section of the process considered and 
V is the velocity of the particle ip, the average been taken over the particle 
distribution. It is clear that the particle is kept in thermal equilibrium as 
far as the term {av) is much bigger that the expansion rate H{t), otherwise 
it decouples from the thermal bath. 

As in the condition of thermal equilibrium the density and pressure of a 
non relativistic species (i.e. T <C m, in appropriate units) is much smaller 
than that of a relativistic one, to a good approximation the two quantities 
take the form: 

^2 



PR = ^9*T\ (1.30) 



and 



7r2 



= Pii/3 = — (1.31) 



where 5* is the total number of relativistic degrees of freedom: 



i=bosons i=fermions 
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gi being the number of relativistic degrees of freedom of each particle specie 
and Ti being its temperature (allowing for deviations from the photon tem- 
perature). 

For temperature T > 300GeV all particles of the Standard Model should 
be relativistic and a value g^: ~ 100 is obtained. For temperatures T ~ 
lOOMeV among the known particles only the electrons, positrons, neutrinos 
(the tree flavour of them) and photons remain relativistic and 5* drops to 
about 10. As the temperature falls below 0.5MeV however also the electrons 
and positrons slow down and a value 5^, ~ 3 is obtained. 

During the radiation-dominated epoch, the use of pR and pR in the Fried- 
mann's equations yield the following useful relations: 

and ^ 

t = 0.301^7^/^ ^^^~(^) sec. (1.34) 

The evolution of the temperature with the scale length R can be under- 
stood in terms of the conservation of entropy S. For this purpose it useful 
to define the entropy density as s = It follows that whenever g'* is 

constant the result T oc R~^ is obtained. The number of effective relativistic 
degrees of freedom for the entropy is defined as: 

g*s = J2 Si (y) + ^ Si , (1-35) 

i=bosons i=fermions 

with notation similar to Eq. 11.321 

In brief the thermal history of the primordial plasma is the following: 
in the first phase the space was filled by an almost homogeneous plasma 
of elementary particles at thermal equilibrium and at very high temper- 
atures (10"^"^ s, more than 10^^ GeV characteristic energy). As the scale 
factor increases the temperature drops and the plasma undergoes several 
phase transitions, most notably the inflationary and GUT (Grand Unified 
Theory, 10~^* s, 10^^ GeV) ones, then the Electroweak phase transition 
(10~^° s, 10^ GeV) and the QCD one (Quantum Chromo Dynamics, 10~^ s, 
10~^ GeV). Among the relics of these phase transitions there are the pri- 
mordial density fluctuations left from the inflationary epoch and the baryon 
content of the Universe after the QCD transition. These eras are followed by 
the nucleosysnthesis era (1 — 200 s, 1 — 0.1 MeV), neutrino decoupling and 
electron-positron annihilation (1 min, 0.5 MeV). Much later the matter and 
radiation have the same density (10^ yrs, 1 eV), and afterwords the electrons 
become bound to the nuclei to form atoms (3 x 10^ yrs, 0.3 eV). Then the 
CMB photons decouple from the plasma traveling freely. From this epoch 
on starts the formation via gravitational collapse of visible structures. 
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The abundance of a particle specie at the decoupling can be estimated by 
properly manipulating Eq ll.291 For this purpose let us define the two vari- 
ables Y = n/s and x = m/T, with m mass of the particle specie considered 
and s the entropy density. Prom the entropy conservation it follows: 

h + 3Hn = sY. (1.36) 

Moreover during the radiation dominated epoch the relation between time 
and temperature obtained above reduces to: 

t = 0.301g:'^'^^^^x'. (1.37) 

Prom Eq ll.291 then it follows: 

dY —(av)s / 9 ON , s 

The exact solution of Eq. 11.381 depends on the cross section a, which in 
turn depends by the particle physics model adopted. However introduction 
the parametrization: 

{av)^a + 6b/x, (1.39) 

valid for non relativistic species, the relic density expressed in terms of the 
critical density assumes a simple form: 

,2 1.07 X 10^ GeV^^ XF 1 

(87rG)V^ .^a + 6b/xF 



3 X 10-27 cm^s-i 
{av) 



fl.41) 



where xp = m/Tp is the x parameter evaluated at freeze-out temperature. 
Por a particle with a given mass, the annihilation cross section has an upper 
bound imposed by the unitarity of the S matrix: {av) ~ 1/m? . This limit 
can be transformed in an upper limit for the DM particle mass by taking 
the DM abundances from the recent WMAP measurements: m < 120 TeV. 
Por more precise estimation of the relic abundances see [22j and references 
therein. 
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Chapter 2 

Dark Matter Particles 



As discussed in the previous chapter the ACDM paradigm needs a compo- 
nent which behaves hke a non-relativistic pressureless dark component. This 
component may well be represented by particle candidates which extend the 
standard model of particle physics at a scale above hundred GeV. It is re- 
markable that such modifications are expected also from a pure theoretical 
reasoning giving good synergy between astrophysical observations and par- 
ticle physics. 

This chapter discusses the fundamental properties of the hypothetical 
new particles giving rise to DM. The possibility of direct or indirect detection 
is also discussed. Specific models are hence presented, highlighting the most 
appealing candidates. 

2.1 Fundamental properties 

Astrophysical and cosmological measurements provide elements that DM 
particles were already present in a non-relativistic state in the early Uni- 
verse. These observations clearly put constraints on the life-time of the 
candidate to be r > 4.3 x 10^^ s. Moreover these particles should interact 
with the already known particles at most weakly, hence excluding charged 
particles (which would not be dark, if not in very specific models excluded 
however by experiments) or particles with color quantum numbers (see [22| 
and references therein for a review). 

It has been proposed that DM may be subject of self interaction. This 
interaction would help in solving the cusp vs core controversy (discussed 
extensively later in this Thesis) for values of the cross section per unit mass 
0-3 ^ a/m < 10^ cm^g~^. The recent observation of the merger of two 
clusters (the so called Bullet cluster) however firmly constrains the cross 
section to a/m < 1 cm^g^^. Other weaker observations further constrain 
the allowed self-interaction, making of it a disfavoured hypothesis. 

Also the self-annihilation has been proposed as a mechanism to reconcile 



17 



18 



CHAPTER 2. DARK MATTER PARTICLES 



the cusp vs core controversy, this mechanism however is excluded by both 
astrophysics (a self-annihilation would produce density cores of same radius 
for different galaxies, in contradiction with observations) and particle physics 
measurements. 

Another important constrain on the DM properties comes from the mod- 
els of stellar evolution. If DM was significantly produced in the interior of 
sun-like stars (due to the high temperature condition), it would change the 
energy loss rate of the sun core modifying all the stellar evolution mechanism. 
This observation provides strong bounds on CDM candidates based on light 
particles such as the axion. Similarly also the BBN measurements provide 
important constraints on the light candidates properties (m < 1 MeV) as 
well as on the decay rates of some heavier particles in some specific models. 

Although severe constraints exist on the DM properties, one of the re- 
quirements of any realistic model is the correct prediction of its abundance 
^DM- The various model can be divided in two classes: the one with thermal 
production (which advocate WIMP candidates) and the other non-thermal 
models (whose prototype is the axion). The thermal models are currently 
more developed and better constrained. They however require some kind of 
weak interaction which can be tested in on going experiments or observa- 
tions. 

2.2 Direct Searches 

The direct searches are focused in detecting the DM direct interaction with 
ground based detectors. The two possible interactions are either with elec- 
trons or with nuclei, of the two however only the second have an acceptable 
sensitivity while the huge background due to natural radioactivity (mainly 
beta decays) makes the first unfeasible. Moreover only DM particles with 
mass sufficiently high can generate a nuclear recoil with detectable energy 
transfer, making of the WIMPs the only acceptable candidates for this kind 
of searches. 

The interaction rate on the detector depends on three quantities: the 
DM flux, the DM-nucleus cross section and the detector mass. Assuming a 
local density of pdm ~ 0.3 GeV cm"^ and a mean velocity of the same order 
of one of the sun around the galactic center {v ~ 220km s^^) the expected 
flux is $ ~ 10'''(GeV/mDM)cm^^s^^, where m,£)M is the DM particle mass 
expressed in GeV. Clearly the big astrophysical uncertainty in the deter- 
mination of the local density affects directly the detection rate prediction. 
Moreover the actual calculation involves not only the mean particle velocity 
but the full distribution of velocities, making the prediction even more un- 
certain. The revolution of the Earth around the Sun however modulates the 
mean velocity according to: 

v{t) = 220 km/s {1.05 + 0.07cos [27r(i - tm)] /I year} , (2.1) 
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where tm is approximately the begin of June. This modulation offers an 
important handle for detecting the DM signal. 

Even bigger uncertainties arise on the cross section side. Not only the 
cross section depends on the the Particle Physics model under study, but also 
large uncertainties arise from the theoretical description of the chosen target. 
In general two kinds of cross sections are studied: the Spin Independent (SI) 
and the Spin Dependent (SD). An important characteristic of the SI coupling 
is that it is coherently enhanced in nuclei according to: 



where A is the atomic number, Mj-ediMN, M^) and Mred(p, M^) denote the 
reduced mass of the WIMP-Nucleus or WIMP-Proton systems respectively. 
The SD coupling normally does not have a similar enhancement, making 
experiments with heavy nuclei far more sensitive to SI interaction in most 
of the cases (although particle physics models with SI coupling suppression 
or SD enhancement exist). 

The detection strategy is based in detecting one or more of the following 
effects generated by DM-nuclei interaction: 

• Ionization: electrons liberated by the atom in primary or secondary 
interactions 

• Scintillation: photons emitted by the de-excitation of excited atoms 

• Heat: phonons generated by the displacement of the nucleus with re- 
spect to the crystalline structure of the detector 

Typically the experiments are sensitive to energies deposited in the detector 
above the keV magnitude. In the last decade several experiments have been 
run, most often with null results. What follows is a brief description of the 
most relevant ones. 

CDMS - Cryogenic Dark Matter Search [23]: this experiment, now 
running deep underground in the Soudan facility, employs crystals of Sili- 
con or Germanium kept at temperatures as low as 10 mK. The detectors, 
known as ZIP detector and featuring the state of the art thin film supercon- 
ducting technology, aims at detecting both the phonons and the ionization 
signals. The combination of the two signals allows a precise constrain on the 
background, especially induced by neutrons. 

Edelweiss - Experience pour DEtecter Les Wimps |^: as for 
the CDMS experiment, the Edelweiss detection technique is based on the 
coincidence of heat and charge detection. In this case however the heat is 
measured by very sensitive thermometric sensor glued on the Germanium 
crystals. 




(2.2) 
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WARP - Wimp ARgon Programme [25]: this experiment, located 
at the Gran Sasso facihty, searches for nuclear recoils in liquid Argon with 
deposited energy in the range 10 — 100 keV by means of both ionization and 
scintillation. The advantage of this technique over the Silicon or Germanium 
detectors relies in the capability of the Argon based detectors to be more 
easily scalable to higher fiductial masses, increasing the sensitivity. 

XENON Dark Matter Programme [26]: as in the WARP case, this 
experiment aims at measuring both the charge and the light signal. Although 
the technical detection details are different, also the XENON experiment has 
the advantage of being relatively easy to scale to high fiductial masses. This 
experiment has recently published one of the most stringent limits on WIMP 
particles. 

CRESST - Cryogenic Rare Event Search with Superconducting 
Thermometers [2^: in this case the detection is based on the combination 
of scintillation and phonon detection. As the active targets are crystals of 
CaW04, the SI interaction is enhanced due to the high mass number of 
tungsten. 

DM-TPC - Dark Matter Time Projection Chamber [2^: this is a 
novel detection scheme based on a low pressure gaseous detector. The exper- 
iment should be able to measure a small track of the recoiled nucleus (which 
should travel few mm), making possible the measurement of the direction 
of arrival of the WIMP particle and hence providing a powerful tool for the 
study of the annual modulation of the signal. 

DAMA [29j : this experiment measures the scintillation in Nal crystals. 
This highly controversial experiment is the only one having reported a sig- 
nal detection. As the experiment is based on only one detection technique, 
the background suppression and the control of other systematic effects are 
more difficult, however the advantage is in the capability of lowering the 
detection energy threshold. Fig. 12.11 shows the annual modulation of the 
detected signal as a function of time for recoil energies between 2 and 4 keV. 
It is interesting to note that the amplitude, period and phase of the modu- 
lation is actually compatible with Eq. 12.11 It is puzzling however that other 
experiments with similar or better sensitivity did not find any signal; it is 
true however that the comparison between different experiments is somehow 
model dependent. The DAMA result, surprising and controversial, will be 
carefully checked by future experiments and certainly the signal detected by 
just one group is not sufficient to claim for a discovery. 

The results of the most sensitive experiments, together with the DAMA 
signal and some theoretical predictions are shown in Fig. 12.21 Clearly most of 
the models predict a cross section several order of magnitude below the cur- 
rent experimental sensitivity, however future improvements in the fiductial 
mass of the experiments together with longer time exposure will hopefully 
improve the situation. The range of masses and cross sections allowed by 
the DAMA signal is model-dependent, however the reconciliation with the 
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Figure 2.1: Annual modulation of the DAMA signal pO]: evidence of direct 
DM Detection? 



other experiments is non-trivial. 



2.3 Indirect Searches 

The Indirect searches of DM are based on astrophysical observations of the 
products of DM self-annihilation or decay. Given the known long lifetime 
of the DM the signal for decay products is suppressed for heavy candidates 
(due to the combination of low number densities and long lifetime) leaving 
only the self-annihilation as most sensitive possible source of a signal. 

In the case of searches via gamma ray observation, the expected flux in 
a detector on Earth is given by: 

dNJ: If f „ . , . o . . . . 

dE^ iirm^ ^ dE^ 2 J^^ Aip Ji ^ .,^ 

where E-y is the photon energy, is the DM particle mass, Atp is the de- 
tector opening angle, {crv)ann is the mean annihilation cross section times 
the relative velocity (of order 10~^^ cm^s~^ for cold WIMP relics from abun- 
dances constraints), Bf indicates the branching fraction in a given channel /, 

is the photon spectrum for a given annihilation channel which depends 
on the DM model and can have both continuum and discrete lines contribu- 
tions, p is the DM density and the integrals are along the line of sight and 
over the detector opening angle. 

The Quadratic dependence on p suggest that the preferred targets for 
indirect searches are the places with higher DM concentrations, like the 
centre of galaxies or galaxy clusters. It must be noticed however that the 
galactic centres are very often sources of strong activities due for example of 
the presence of a Black Holes or other compact objects enhancing the over- 
all background. Moreover the large uncertainty on the DM density reflects 
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Figure 2.2: Search of SI DM interaction: limits, DAMA signal and theoretical 
predictions. 



directly on the flux predictions making the searches extremely difficult (al- 
though possible enhancements due to local DM over-densities are possible). 

Another possibility is to pursue indirect search by looking at charged 
particles such as positrons or antiprotons, in this case however the galactic 
magnetic is such that the direction of arrival of the particle does not reflect 
the production point and the only observable if an excess of antimatter with 
respect to the expected background due to ordinary cosmic rays (which also 
suffer from big uncertainties). 
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What follows is a brief partial review of the most important facilities 
looking for indirect signals: 

XMM Newton and Chandra [3H[32]: these are two satellites operated 
by the European Space Agency (ESA) and by the National Aeronautics and 
Space Administration (NASA) respectively. They have both imaging and 
spectroscopic capabilities in a photon energy window between 0.1 keV and 
10 keV approximately. The searched signal is a narrow line not explicable 
in terms of weak known physics processes and originating either from DM 
decay or self annihilation. Their observations put important limits on the 
medium mass DM candidates (such as the sterile neutrinos). 

Integral fSSj: this ESA observatory operating a window of energy of 
gamma rays between 15 keV and 10 MeV approximately is complemented 
by optical instrumentation. This observatory may detect a signal of DM 
as a new narrow line, as an excess of 511 keV photons due to positrons 
annihilation. The mission has actually published a claim of possible DM 
detection discussed below. 

Compton Gamma Ray Observatory - CGRO [3^: this observa- 
tory, together with Hubble and Spitzer, is one of the most important re- 
search projects of NASA. Two instruments on board made important DM 
searches: COMPTEL (Imaging Compton Telescope) operating in an energy 
range of 0.75 MeV and 30 MeV and EGRET (Energetic Gamma Ray Exper- 
iment Telescope) operating in the window 20 MeV - 30 GeV. The EGRET 
telescope in particular provided important limits both on the DM properties 
and possible signals of detection discussed below. 

AGILE and GLAST [351 EB]: these are the two recently lunched 
gamma ray observatories. The first operates in the energy windows 30 MeV— 
50 GeV and 10 — 40 keV, allowing both gamma and X-ray measurements. 
The second one has full coverage of the window 10 keV — 300 GeV with both 
large opening angle and excellent sensitivity. GLAST, with its unprecedented 
sensitivity, is certainly the best observatory to look for indirect searches of 
WIMP annihilations. 

CANGAROO, HESS, MAGIC and VERITAS [37l [381 [39l [40]: 
these are ground based facilities observing Ultra High Energy gamma rays 
with energies above 100 GeV approximately. The detection technique is 
based on the measurements of the Cherenkov light emitted by electromag- 
netic showers in the upper atmosphere. Constraints (not very rigid yet) on 
self annihilating WIMP particles in nearby halos have been provided. 

AMANDA, ICECUBE and ANTARES dH ill HS] : these are High 
Energy Neutrino observatories. The indirect search of DM with this detector 
is based on the assumption that a high density of WIMP particles would 
accumulate at the core of the Sun or of the Earth due to a combination of 
the elastic scattering of the particles with the Sun or Earth material followed 
by a gravitational capture. The local high density of DM would enhance 
the self-annihilation which may proceed through a channel which includes 
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neutrinos in the final state. This neutrinos would then easily escape from 
the Sun or Earth core allowing the detection on the earth surface. No signal 
has been found so far. 

PAMELA and AMS [111115]: these are two satellites which aims at 
measuring the spectra and properties of primary cosmic rays in the GeV 
region. The presence of an unexplained excess of antimatter (either antipro- 
tons or positrons) in the primary cosmic rays can be interpreted as a signal 
of self-annihilating DM. It has to be mentioned that few years ago the HEAT 
experiment found evidence of an excess of positrons with energies of about 
7 GeV and the signal has been confirmed later by AMS-I. The interpretation 
of the signal as coming from DM annihilations however is problematic due to 
the much lower excess predicted by the WIMP models. However new data 
from PAMELA should clear the uncertainty soon. 

Up today several claims of indirect DM detection has been made, some- 
times in conflict with each other or with other measurements. The most sig- 
nificant are however: the positron excess measured by HEAT, the 511 keV 
fine excess measured by INTEGRAL, the EGRET Diffuse Galactic Spec- 
trum, the EGRET Diffuse Extragalactic Spectrum and the so called WMAP 
Haze (an excess of microwave emission around the center of the Milky Way). 

The INTEGRAL signal is many order of magnitude above the expected 
signal from secondary positrons due to cosmic rays and is approximately 
spherically symmetric with a full width half maximum of about 6°. Astro- 
physical interpretations of the signal are difficult and several interpretations 
due to indirect DM detection have been proposed. 

The EGRET Galactic and Extragalactic Spectra are shown in Fig. 12.31 
The Galactic measurements show an excess of photons with energies in the 
range 1 — 10 GeV approximately. The interpretation of this signal as DM de- 
tection however is questionable due to the mismatch between the knowledge 
of the halo density distribution and the directional variation of the signal, 
moreover a large amount of secondary antiprotons would be expected in con- 
trast with observations. The Extragalactic excess measured by EGRET with 
energies above 10 GeV can also be interpreted as a DM signal. In order for 
this interpretation to be valid however the DM halos have to be very cuspy 
for most of the galaxies BUT far less cusped for the Milky Way, a rather odd 
situation indeed (beside being in confiict with other observations). 

The last of the above mentioned claims of indirect detection, the WMAP 
Haze, has been proposed to arise from synchrotron radiation emitted by 
relativistic positrons or electrons generated by DM annihilations. If this is 
correct however an associated prompt gamma ray emission should be within 
of the recently lunched GLAST experiment. 

To conclude this section a final remark is necessary: several claims have 
already been proposed as indirect DM detection, all of them however com- 
peting with other Astrophysical explanations. It seems clear then that a con- 
vincing evidence of indirect detection must come from several complementary 



2.4. THE ZOO OF CANDIDATES 



25 




Figure 2.3: On the left: EGRET measurement of the Diffuse Galactic Spec- 
trum on a portion of the sky as evidence of DM (see ^6] and references 
therein). On the right: EGRET measurement of the Diffuse Extragalac- 
tic Spectrum as evidence of DM (with different characteristics from the left 
panel) . 

measurements (at different energies or with different particles) although all 
the claims have to be properly studied and possibly eliminated. 

2.4 The zoo of candidates 

Since the early years following the discovery of DM halos surrounding galax- 
ies particle physicists tried to develop models which included a DM can- 
didate. Today almost every proposed extension of the Standard Model of 
Particle Physics (SM in the following) in a way or another includes a can- 
didate. Most often these candidates are WIMPs produced thermally in the 
early Universe, and in this case the long lifetime is ensured by including in 
the model a symmetry which forbids the DM decay. Other models with non- 
thermal production however also play an important role in the discussion of 
the extensions of the SM. 

What follows is a brief discussion of the two most attractive candidates: 
the supersymmetric models and axions. The list of other proposed candidates 
is however very long and includes: sterile neutrinos, minimal DM models. 
Little Higgs models, Kaluza-Klein particles, wimpzillas, CHArged Massive 
Particles (CHAMPs), brane- world DM and many others. 

Supersymmetric candidates: In the SM bosons and fermions play 
two different roles: the bosons act as mediators of fundamental interactions 
while the fermions are the elementary constituents of matter. It is natural 
therefore to ask whether a symmetry exists between these two particle classes 
providing a sort of unified picture. This boson-fermion symmetry is called 
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super Symmetry (SUSY). 

While a review of the SUSY theory is beyond the scope of this Thesis, 
the discussion here will be concentrated on the most common SUSY DM 
candidate (the neutralino), referring the reader to some excellent reviews 
and didactic materials available in the literature (see e.g. ^7] and references 
therein). 

The benefits of the SUSY models include not only a suitable DM candi- 
dates, but help in solving the so called "hierarchy problem" (the difference 
between the electroweak and the Planck scales) and provide a mechanism 
for the unification of the gauge coupling of the SM at a Grand Unification 
Scale. 

In the Minimal Supersymmetric Standard Model (MSSM) each boson 
(gluons, and B) is associated with a fermion (the gluinos, winos and 
binos), the quarks and leptons are associated to scalars called squarks and 
sleptons and the Higgs sector is composed by two Higgs doublets associ- 
ated with spin 1/2 higgsinos. Another ingredient of the MSSM is a discrete 
symmetry called i?-parity where to each particle is associated a conserved 
quantum number defined as i? = (_i)3^+^+2s^ where B and L are the bary- 
onic and leptonic number respectively, while s is the particle spin. Clearly 
all the SM particles have R = 1 while the SUSY partners have R = —1 and 
the i?-parity conservation implies the stability of the lightest SUSY partner. 

A consequence of SUSY is that the mass of each particle must be equal 
to the mass of its super-partner, otherwise SUSY is broken. Clearly the mass 
degeneracy predicted by SUSY is not observed in Nature and several SUSY- 
breaking mechanisms have been proposed. It is evident that the original 
SUSY idea led to the formulation of a model with some attractive feature 
(like grand-unification of the couplings) at the price of an enormous increase 
of the number of free parameters. It is common however to try to reduce the 
number of free parameters by introducing some kind of additional condition 
like the unification of the gaugino masses at GUT scales or some universality 
of the couplings. The most common resulting models are the Constrained 
MSSM (CMSSM) or the minimal Super Gravity (mSUGRA). 

As the MSSM (with its variants) received a lot of attention in recent 
years both from the theoretical and experimental community important con- 
straints exist on the parameter space. Fig. 12.41 shows for the CMSSM in the 
parameter space (mo; rni/2) (universal sfermion mass in the vertical axis and 
universal gaugino mass in the horizontal one) for two definite values of the 
other CMSSM parameters the region allowed by the WMAP measurements 
(turquoise) together with bounds from accelerator measurements or other 
observations. An immediate observation which can be drawn from these di- 
agrams is that the parameter space compatible with the cosmological bounds 
is well constrained and an important portion of this region is not compatible 
with the other measurements. Clearly all the bounds can be relaxed in a less 
constrained MSSM, in this case however the huge number of free parameters 
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Figure 2.4: On the left: the ("1-1/2, "t-o) planes for tan/? = 10 and fi > 0, 

assuming ^0 = 0,mt = 175 GeV and mi){mh)^^ = 4.25 GeV. The near- 
vertical (red) dot-dashed lines are the contours nih = 114 GeV, and the near- 
vertical (black) dashed line is the contour m^± = 104 GeV. Also shown by the 
dot-dashed curve in the lower left is the corner excluded by the LEP bound 
of mg > 99 GeV. The medium (dark green) shaded region is excluded by 
6 — > 57, and the light (turquoise) shaded area is the cosmologically preferred 
region. In the dark (brick red) shaded region, the LSP is the charged fi. The 
region allowed by the E821 measurement of at the 2-a level, is shaded 
(pink) and bounded by solid black lines, with dashed lines indicating the 1-a 
ranges. On the right: tan/3 = 50 [47| . 

makes the model less attractive. 

Axion: One of the open problems of the SM is that the gauge theory 
responsible for the strong interaction foresees the possibility of a strong CP 
violation (see |i48j and references therein). This strong CP violation however 
is not observed in Nature, hence the so called strong CP problem arises. One 
of the possibility for its solution is that Nature respects a symmetry, called 
Peccei-Quinn symmetry, which allows the restoration of the CP conservation 
in the strong sector. The PQ symmetry however should be spontaneously 
broken giving rise to a new Nambu-Goldstone called Axion. 

The Axion has a specific property of being coupled to two photons as: 

Ca^ = -ga-y^-^a, (2.4) 

where E and B are the electric and magnetic fields respectively, a is the 
axion field and the coupling constant (7a7 is related to more fundamental 
parameters of the theory such as the axion mass: Qa-y oc rua. 

The coupling plays a fundamental role in the Axion searches. Axions in 
fact can transform into photons when propagating in an external magnetic 
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field in a way similar to neutrino oscillations. As visible in Fig. 12.51 both 
astrophysical and laboratory measurements impose strong constraints on the 
properties of a hypothetical axion with mass as low as about 1 eV. Prom 
the current DM search point of view the on-going Axion DM eXperiment 
(ADMX) is searching for a signal in a mass region close to the ;^eV scale. An 
axion with a mass around this scale in fact is a theoretically well motivated 
cold DM candidate. 

While axions with a mass above the eV scale would be produced thermally 
in the early universe (and hence would be a hot DM candidate similar to 
the neutrinos) for masses lower than the /ueV the production would be non- 
thermal and linked to the Peccei-Quinn phase transition by the so-called 
misalignment mechanism. The relic axion density can then be calculated 
according to: 



where — vr < 0i < vr is the initial "misalignment angle" relative to the CP- 
conserving position and fa oc is the Peccei-Quinn scale. 

In this case an axion with rua ~ lO^eV would provide a CDM density 
in agreement with WMAP measurements, however, this number sets only a 
crude scale of the expected mass for axion DM, with uncertainties coming 
both from the particle physics and the cosmological models. It has to be 
mentioned that in the non thermal axion production mechanism the effective 



temperature! 

2.5 Concluding remarks 

The dedicated searches of the DM particle candidate have seen an important 
boost in recent years with relevant and costly experiments been planned and 
executed. The success of these searches however crucially depends on our 
capability of predicting the signal expected for each particle model. This 
capability in turn relies on our knowledge of the DM distribution in bound 
objects. The mass distribution in galaxies is then the theme of the next 
chapter. 




(2.5) 




an extraordinary low 
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Figure 2.5: 



Axion limits and foreseen search ranges [48j . 
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Chapter 3 



The mass distribution in Spiral 
Galaxies 

The presence of large amounts of unseen matter in spiral galaxies with a 
distribution different from that of stars and gas is well established. The pri- 
mary observational evidence for the existence of DM, under the assumption 
of Newtonian gravity, comes from optical and 21 cm RCs of spirals which do 
not show the expected Keplerian drop-off at large radii but remain increas- 
ing, flat or gently decreasing over their entire observed range [49l [501 ISTl [52] . 
The invisible mass component becomes progressively more abundant at outer 
radii and for the less luminous galaxies [531154] • The distribution of matter in 
disk systems has become a benchmark for the present understanding of the 
process of galaxy formation and, with the help of the available observational 
tools, crucial questions can be addressed: 

• has the dark matter an universal distribution reflecting its very Nature? 

• how and why the dark-to-luminous mass ratio and other physical quan- 
tities vary in objects of different Hubble type? 

• how dark matter affects the fate of the universe? 

It is well known that numerical simulations performed in the ACDM sce- 
nario predict a well-defined density profile for the virialized halos surround- 
ing and hosting the galaxies. This profile leads to structural properties of 
galaxies [3] that are in strong disagreement with observations. Moreover the 
mechanism of galaxy formation, as currently understood, involves the cool- 
ing and the condensation of HI gas inside the gravitational potential well 
of DM halos. Part of the condensed gas then transforms into stars which 
reheat the former by the feedback of SN explosions. It is clear then that a 
mistake in the model of the halo potential has a deep impact on the complex 
dynamics of stars and gas. 
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Figure 3.1: On the left orbital velocity of the planets of the solar system as 
a function from the distance from the Sun (Copyright ©Addison Wesley). 
This picture shows the typical Keplerian falloff planetary systems as well as 
demonstrating that on these scales the gravitational potential is dominated 
by the Sun mass with negligible contribution from DM. On the right circular 
velocity of the NGC 1090 disk. It clearly shows the absence of any Keplerian 
falloff and hence revealing the presence of DM. 



It is widely accepted that the mass distribution of spiral galaxies, which 
can be derived from observations, bears the imprint of the Nature and the 
cosmological history of DM. Moreover it reflects the interaction (possibly not 
only gravitational) between dark and luminous matter. 

Although the DM presence is widely accepted a strong debate about its 
spatial distribution is ongoing and very little is known about its Nature. 

After a brief historical introduction, RCs are described as main tracers of 
gravitational field in galaxies. Then their most important properties, includ- 
ing their slope as well as scaling relations, are presented. These properties 
are well described by the Universal Rotation Curve (URC) paradigm. 

The predictions of the halo properties from numerical simulations is the 
subject of a later paragraph, which is followed by an extensive discussion of 
the the comparison with the observational properties of the inner distribution 
of DM. 



3.1 Historical Introduction 

The use of galaxy kinematic as a tool for studying the mass distribution has a 
long history (see ^5] and references therein) starting in the second decade of 
the XX century by works of Slipher [56| and Wolf |57j . In particular it is due 
to Slipher the discovery that the Andromeda (M 31) galaxy is approaching 
the Milky Way with a speed of order 100 Km/s and its disk is rotating around 
its center with a steeply rising velocity in the inner region. Only more than 20 
years later however the measurements of Babcock ^8] and Oort [59] where 
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precise enough to show that the total mass distribution in spirals is not 
simply proportional to the distribution of light emitted by stars as it was 
expected. It has to be mentioned that in 1933 Zwicky, in a failed attempt of 
analyzing the coma cluster dynamics, made the hypothesis of the existence 
of a mysterious dark component. He was not aware however at that time of 
the importance of its fortuitous statement [60]. The technological advances 
after the second world war allowed astronomers to routinely detect the RCs 
of tens of galaxies. In 1972 Whitehurst &z Roberts [6l] found an anomalous 
high velocity of neutral hydrogen gas around M 31 giving a first hint of non- 
keplerian fall at large radii. Only with the work of Rubin [51] however it was 
clear that "the conclusion is inescapable that non-luminous matter exists 
beyond the optical galaxy". This work opened the so called "dark matter 
problem" in galaxies. The work of Rubin, based on optical observations, 
received an important confirmation and extension by the work of Bosma 
[5U] with HI measurements up to larger galactocentric radii. Clearly at 
that time the Standard Model of Particle Physics was not complete yet so 
the nature of DM component was an open question. However neutrinos 
provided a viable candidate (but soon clearly understood to be excluded). 
The phenomenological analysis of the RCs had an important step further in 
1988 when Persic & Salucci ^53j found a general trend with larger baryonic to 
dark mass discrepancy in fainter galaxies and vice- versa. By that time it was 
clear that the disks of spiral galaxies are embedded in a much bigger spherical 
halo whose nature still remained unclear. Under the assumption of DM being 
a cold collisionless particle in 1996 Navarro, Frenk & White [3] developed 
a computational model for the formation of the halo and gave a simple 
parametrization of the halo mass distribution. On the phenomenological 
side instead the study of more than a thousand of galactic RCs by Persic, 
Salucci & Sersic [62] revealed that they can be well-represented by a URC, 
function of the galaxy luminosity. 



3.2 Rotation Curves as gravitational field tracers 

A RC of a spiral galaxy can be defined as the diagram of the circular velocity 
as a function of the galactrocentric distance, and is the fundamental probe 
of the behaviour of the gravitational potential of the system. 

The mass distribution in a spiral can be modeled as the sum of three 
discrete components: a halo of DM, a disk of stars and gas, and a stellar 
bulge. The halo and the bulge are assumed to have a spherical distribution 
while the disk is approximately as infinitesimally thin, with the centre of the 
three distributions being coincident. The total gravitational potential (ptot 
can then be decomposed as 

4>tot = 4>DM + 4>disk,stars + 4>disk,gas + 4>bulge- (3.1) 

Assuming that the components of the disk have circular orbits with velocity 
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Figure 3.2: Illustration of the measurement of a RC. On the left: points on 
the major axis equidistant from the centre are red-shifted or blue-shifted by 
the same amount. On the right: the RC is obtained from the Doppler shift 
measurements along the major axis. 

ytotii") at a radius r the following relations hold: 

^ot('^) = ~ ^DM + ^Jsfc, stars + ^Jsfc,sas + ^bulge^ (3-2) 

where we defined V^j^j = r d(j)DM/dr, and similarly for the disk and bulge. 

RCs are obtained by measuring the Doppler shift of absorption or emis- 
sion lines of disk material (see Fig. 13. 2p . The most commonly used tracers 
are Ha emission lines obtained by optical spectroscopy, that gives the kine- 
matics of the inner part of the galaxy (stars), and neutral hydrogen HI 
("21-cm line") obtained by radio measurements, that extends up to a larger 
radii. Radio observations have an angular resolution bigger then optical, but 
better spectral resolution corresponding to smaller errors in the velocity. 

It is possible to estimate the contribution of disk and bulge to the total 
gravitational potential from the measurements of their mass surface densities. 
These in turn yield the contributions Vj-^i^^^^^^^, V^-^^ ^^^ and V^^^g^. 

The stellar mass distribution is given by its luminosity distribution mul- 
tiplied by a mass-to-light ratio, which is assumed to be constant within each 
bulge/disk component. Note that from optical measurements it is difficult to 
disentangle the mass surface density of the disk and of the bulge in the inner 
region of the galaxy. For this reason the best mass models are obtained in 
galaxies with negligible bulges. Moreover more complex mass models have 
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been tried introducing more components. However these additional compo- 
nents increase the degeneracy between free parameters without adding much 
physical information. 

The situation is different for HI measurements where the surface lumi- 
nosity density distribution T,gas gives a direct measurement of the gas mass. 
The halo mass distribution can be either parametrized by a theoretical or 
an empirical model, or derived from the observed RC inverting Eq. 13.21 and 
using appropriate models for bulge and disk. 

Optical observations show that very often the stars in the disk follow the 
exponential Freeman profile |63| 

where Md is the disk mass and Rd is the scale length, the latter being 
measured directly from the observations. It is useful to define the optical 
radius (the radius enclosing 83% of the total light, see [62]), Ropt = 3.2 
as the "size" of the stellar disk. In the same way the stars in the bulge 
very often are distributed according to the Sersic mass density profile (e.g. 
[M] and references therein), which yield the following surface mass density 
profile: 

where Mh is the total projected mass, r is the projected spherical radius. Re 
is the effective radius, n is the index of the profile, a ~ 2n— 0.324 and r[2n] is 
the complete gamma function. The index n is associated with the curvature 
and the concentration of the profile; n = 1 corresponds to an exponential 
profile, while the classical de Vaucouleurs profile is obtained for n = 4. 

From the Poisson equation and using cylindrical coordinates, the poten- 
tial due to disk material reads 

poo p2iT jn 

^disk{r) = -G dr' r'T.disk{r') , (3.5) 

5^disfc('^') is the surface density distribution of the stars in the disk $]£)(r'), 
given by (13. 3|) . or of the gas Y^gasir'), given by an interpolation of the HI data 
points up to the last measured point. Having the stars a simple distribution, 
equation (|3.3p can be integrated in terms of Bessel functions and results in 
the usual expression [63| : 

r2 , . GMd 2r.fX 



Vdisk,stars{r) = -^^x B (^-j , (3.6) 



where x = r/R^, G is the gravitational constant and the quantity B = 
IqKo — hKi is a combination of Bessel functions. 



36 CHAPTER 3. THE MASS DISTRIBUTION IN SPIRAL GALAXIES 



-24 
-23 
-22 
-21 
S"-20 
-19 
-18 
-17 
-16 - 



1.0 



opt 



opt 



.•f •..<^•.^f•^■ 
j'.>./i • 



1.6 



2.0 



2.2 



logV 



2.4 



2.6 



Figure 3.3: The Radial TF relation [67] . Each one of the six relations is 
indicated with different colours. 



For a spherically symmetric bulge distribution, one has simply V^{r) = 
Gm{r)/r, where m(r) is the mass interior to radius r. Following ^5], the 
bulge mass gives 

m{r)= dr'2Trr'T,B{r')+ dr [sin-^{r/r') -r{r''^ -r^)~'^^^]4r' T.B{r'). 

Jo Jr 

(3.7) 

The above mass model of RCs is valid under the hypothesis of circular 
motions. The issue of testing this hypothesis then arises. Tully & Fisher 
|66j discovered that the maximal rotational velocity Vmax of a spiral galaxy, 
measured by the full width at half-maximum of the neutral hydrogen 21- 
cm line, correlates with the galaxy luminosity by means of a power law 
of exponent a ~ 4. This equivalently reads M = a logVmax + ^, where 
M is the absolute magnitude in a certain band and 6 is a constant. This 
relation is a powerful tool to determine the distances of galaxies and to 
study their dynamics |53j. The rotational velocity reflects the equilibrium 
configuration of the underlying galaxy gravitational potential. In a recent 
work it has been found a new Tully Fisher relation for spirals holding at 
different galactocentric radii, called Radial Tully-Fisher relation [67] : 

Mband = an log Vn + bn, (3.8) 

where Vn = Vrot{Rn), and a„, 6„ are the slope and zero-point of the relations, 
with Rn = {n/5)Ropt- This relation proves that the rotation velocity of 
spirals is a good measure of their gravitational potential (see Fig. 13.31) . 
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More specifically, the fact that in any object and at any radius, the rotation 
velocity can be predicted just by the galaxy luminosity implies that non 
circular motions are generically negligible. 
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Figure 3.4: RCs of spiral galaxies of different luminosities. 



3.3 Rotation Curves are not flat! 

The information about the distribution of luminous and DM in disk galaxies, 
as well as on the correlations among the main parameters that characterize 
both components, offers a fundamental clue to understand how galaxies form 
and evolve, what role DM plays in these processes, and what imprints DM 
leaves about its nature. It is crucial to remark the observational fact that 
the RCs are not asymptotically flat (see a representative sample of RCs in 
Fig. 13.41) . as it is assumed in a huge number of papers. When in the late 
1970s the phenomenon of DM was discovered [50l EI] a few truly flat RCs 
were highlighted in order to rule out the claim that non Keplerian velocity 
profiles originate from a faint baryonic component distributed at large radii. 
At that time a large part of the evidence for DM was provided by extended, 
low-resolution HI RCs of very luminous spirals (e.g. [50]) whose velocity 
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Figure 3.5: Logarithmic gradient of the circular velocity V ys B absolute 
magnitude and vs log V{Ropt) [7T| . 

profile did show small radial variations. 

The increase in the quality of the RCs soon leads to the conclusion that 
baryonic (dark) matter was not a plausible candidate for the cosmological 
DM and that the RCs did show variation with radius, even at large radii. 
Later numerical simulations in the CDM scenario also predicted asymptot- 
ically declining RCs [3]. The flat RC paradigm was hence dismissed in the 
90's (e.g. [531 EHl EH]). Today, the structure of the DM halos and their 
rotation speeds is thought to have a central role in Cosmology and a strong 
link to Elementary Particles via the Nature of their constituents (e.g. |70|). 
and a careful interpretation of the spiral RCs is considered crucial. 

It must be noticed that the circular velocity due to a Freeman stellar disk 
has a flatfish profile between 2 and 3 disk scale-lengths which implies that 
a flat RC is not necessarily a proof for the existence of DM. Its most solid 
evidence instead originates from the fact that even in very faint galaxies the 
RCs are often steeply rising already in their optical regions. 

A quantitative analysis on the issue is shown in [7T], where the concept of 
RC logarithmic slope, defined as V = {dlog V/dlog R), is used. By plotting 
the logarithmically slope at the optical radius for a huge sample of galaxies 
|62| . see Fig. 13.51 it is clear that: 

-0.2 < V < 1, 



i.e. it covers most of the range that a circular velocity slope could take [-0.5 
(Keplerian), 1 (solid body)]. Notice that a flat RC means V = 0, while 
in the case of no DM the self-gravitating Freeman disk lead to V = —0.27 
at 3 Rd- It is also important to notice the strong correlation between the 
rotation shape (V) and the galaxy luminosity [62j[72l[73] (see Fig. 13. 5|) . 
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The incredibly amount of theories that either imply or assume the exis- 
tence of an observational scenario in which the RCs of spirals are asymptot- 
ically flat, is clearly in contradiction with observational evidences. 

3.4 The Universal Rotation Curve 

The studies of spirals of type Sb-Im in the '90, pioneered by [H] and fur- 
ther developed by [H2], led to the remarkable observation that the RCs of 
these objects present universal properties well correlated with other galactic 
properties like the disk mass or the virial mass. These works led to the con- 
struction of the so called "Universal Rotation Curve" VuRciR, P) [75], i-e. a 
function of the galactocentric distance R tuned by the chosen parameter P 
(e.g. the virial mass). Three different coordinate systems are normally used 
to measure the radius: the physical coordinate R, the radius expressed in 
terms of the scale length R/Rd or in terms of the DM characteristic length 




50 - : 
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Figure 3.6: The URC with the radial coordinate in physical units [73]. Each 
curve corresponds to Mmr = lO^^lC^/^M©, with n = 1 . . . 9 from the lowest 
to the highest curve. 

In order to develop the URC, 11 synthetic curves Vcoadd{R/ Ropti^i) 
were built by selecting 616 RCs of galaxies with negligible bulge or HI disk, 
subdividing them in 11 groups spanning in total the /-band luminosity range 
— 16.3 < M/ < —23.4. Each luminosity bin then contain about 1500 velocity 
measurements (from different galaxies by with similar luminous properties) 
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Figure 3.7: The URC, normalized at its virial value as a function of nor- 
malized radius R/Ryir [75J. Each curve, from the highest to the lowest, 
corresponds to M^r as in Fig l3.6l The bold line is the NFW velocity profile. 

which are arranged in radial bins of size 0.3 Rd up to ~ 4i?£) and then 
coadded. The synthetic curves obtained are then free from most of the 
observational errors and non-axisymmetric disturbances present in individual 
RCs, smooth and with a very small intrinsic variance. The properties of these 
curves are then found to strongly correlate with luminosity (see also fT6\). 
The additional data used in the URC are the empirical relationship between 
RC slope at 2 Ropt and log Vopt (see [62]) and the halo virial velocity Vmr = 
{GMmr / RvirY^'^ , obtained from the disk mass vs virial mass relationship 



The URC paradigm, which states that the halo or disk mass determines at 
any radii the circular velocity of any spiral by means of the URC function, 
is the observational counterpart of the NFW velocity profile obtained by 
numerical simulations. 

The URC function is modeled as the sum in quadrature of two terms: 
^URC ~ ^URCD + ^URCH i where V^^Qj^ represent the disk contribution and 
^URCH the DM halo. 

The disk contribution is given by Eq. 13.61 For the DM term it is assumed 
the empirically Burkert profile [^, a cored distribution that can converges 
to a NFW one at outer radii: 



m 



p{R) 



(3.9) 



{R + ro) (i?2 + 
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Figure 3.8: Best disk-halo fits to the URC (dotted/dashed line: disc/halo) 
[751. 



ro is the core radius and po the central density density. Then: 

yi.cu(m ^ 0.4 c Mt{,„(i + ii) - t.n-. (^) . i,„[i . (^)^] } . 

(3.10) 

The URC function then has three free parameters po, '^o, that are 
obtained from fitting Vcoadd and the other data specified above. In Fig. 13.61 
the URC function Vij^c(R\Mmr) is shown expressing the radius in physical 
units and identifying the objects by the halo virial mass. Each line refers to a 
given halo mass in the range 10^^ Mq — 10^^ Mq. The halo mass determines 
both the amplitude and the shape of the curve. Note however that the 
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contribution of the baryonic component is negligible for small masses but 
becomes increasingly important in larger structures. In Fig. 13.71 the URC 
VuRciR/ Rvir] Mmr) is shown as a function of the radial dark coordinate 

1 /3 

R/Rmr and is normalized by Vmr oc Mj^^ . 

The URC shows that the DM halos and stellar disks are both self-similar, 
but the whole system is not, likely due to the baryons collapse that have 
broken it in the innermost 30% of the halo size. 

RCs are critically not flat: their RC slopes take all of sort of values from 
that of a solid-body system (i.e. +1) to that of an almost Newtonian point- 
mass (i.e. - 1/2). The maximum of the RC occurs at very different radii, for 
galaxies of different mass, viz. at ~ 2R£) for the most massive objects and 
at ~ IQRd for the least massive ones. 

The existence of systematical properties of the mass distribution in spi- 
rals was first claimed by [53] and then successively confirmed by independent 
works [Ml ESI EH EH]. In order to understand the whole process of cosmo- 
logical galaxy formation we must take into account the rich scenario of the 
dark-luminous interplay occurred in galaxies. 

In detail, the mass distribution in Spirals, as carefully obtained in [75], is 
obtained by mass modeling two very different and complementary kinemat- 
ical set of data a large number of individual RCs of objects of different 
luminosity and h) the URC (see Fig. 13.8(1 . The noticeably very similar re- 
sults obtained from these two different sets of data strongly indicates their 
robustness and reliableness. A clear scenario of the mass distribution then 
emerges (see Fig. 13. 9p : 

• The stellar disk dominates the galaxy's inner region out to the radius 
at which the DM halo contribution starts to take over the stellar one. 
This sets the properties of the Radial Tully Fisher relation and yields to 
the paradigm of the Inner Baryon Dominance: the inner observed RC 
that can be accounted by the stellar matter alone are indeed saturated 
by this component. 

• At any radii, galaxies with lower luminosities have progressively more 
proportion of DM i.e. a larger dark-to-stellar mass ratio. In detail, 
the disk mass is oc M^^^ at small halo virial masses (e.g. M^r = 
lO^^M©) and oc M^ir at larger masses (e.g. M^ir = W^Mq). The 
baryonic fraction is always much smaller than the cosmological value 
Qi^/Qmatter — 1/6, and it ranges between 7 x 10^^ to 5 x 10~^ in line 
with is the well-known evidence that SN explosions have removed (or 
made never condense) a very large fraction of the original HI material. 

• Smaller spirals are denser, with the central density spanning 2 order of 
magnitudes over the mass sequence of spirals. 

• The structural parameters of the mass distribution, po, Md, M^, vq 
are remarkably all related, see Figs. 4 and 11 of |62| . 
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Figure 3.9: A summary of the empirical scaling relations between the struc- 
tural parameters of the mass distribution [75]; top left: stellar disk - hallo 
mass; top right: disk mass - scale-length; bottom left: disk mass - halo central 
density; bottom right: halo mass - core radius. 



• The stellar mass-to-light ratio is found to lie between 0.5 and 4. The 
values of disk masses derived as above agree very well with those ob- 
tained by fitting their SED with spectro-photometric models |80| . 

• The HI component is almost always below the kinematical detectably. 
However, in low mass systems it cannot be neglected in the baryonic 
budget since it is more prominent than the stellar disk. 



3.5 Dark halos from simulations 

In the standard picture of galaxy formation, DM halos provide the frame- 
work for the formation of luminous galaxies (e.g., [HH [82l [83]). The DM 
halos are assumed to form hierarchically bottom-up via gravitational am- 
plification of initial density fluctuations. The halos carry with them gas, 
which eventually cools and contracts to form luminous disk galaxies at the 
halo centres. The halo profile has a direct dynamical role in determining the 
observable RC of the disc. It also affects gas cooling and in-fall and there- 
fore the structural properties of the resultant disc, such as size, luminosity 
and surface brightness. In the 1970s numerical simulations were developed 
and used to understand the mechanisms of gravitational clustering, and the 
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evolving quality of the codes, together with the increasing resolution and 
computational power, made them the preferred tool to study the formation 
of Cold DM halos. The success of numerical simulations in reproducing the 
observed dynamical properties of galaxies and larger systems depends on 
the scale investigated, and there is no agreement about the actual shape of 
DM halos and the mass distribution of substructures, due to inconsistencies 
between the results of simulations and observations; however, simulations in- 
deed reproduce well the mechanism of hierarchical clustering, and the latter 
enjoys a much broader consensus in being the actual process responsible for 
structure formation. 



1 2 3 4 1 2 3 4 




-10 1 2-10 1 2-10 1 2 



Log r'adius 

Figure 3.10: Density profiles of simulated halos in different cosmologies [84]. 
In each panel, the lower-mass halo is represented by the leftmost curve; 
the solid smooth curve is the NFW fit. Left panels: Standard CDM model 
(A = 0). Right panels: ACDM model. In each panel the varying cosmological 
parameters are specified. Radii are in kiloparsecs {scale at top); the arrows 
indicate the softening length in each simulation. 

The most evident property of halos born through hierarchical clustering 
is the self- similarity: no matter the mass scale, they all belong to a one- 
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parameter family of curves, known as Navarro, Prenk & White [3j[84] profile 

PNFw{r) = — 2' (3-11) 

{r/vs) (1 + r/rs) 

where is a characteristic inner radius, and ps the corresponding inner 
density. The outer, virial radius Ryir, of a halo of virial mass Myir, is defined 
as the radius within which the mean density is A^j^ times the mean universal 
density pu at that redshift: 

47r 3 

The associated virial velocity is defined by V^^^ = GMmr/Rvir- The one- 
to-one relations between the three virial parameters are fully determined by 
the background cosmology. The virial over-density Avir is provided by the 
dissipationless spherical top-hat collapse model [85j[86]; it is a function of the 
cosmological model, and it may vary with time. For the Einstein-deSitter 
cosmology, the familiar value is A^r — 178 at all times. For the family of 
flat cosmologies {Qm + ^^A = 1), the value of Ayir can be approximated by 
[87] A^ir ~ (187r2 + 82x - 39x'^)/n{z), where x = n{z) - 1, and n{z) is the 
ratio of mean matter density to critical density at redshift z. In the ACDM 
cosmological model (Qm = 0.27), the value is Aj;ir{z = 0) ~ 360. 

An associated useful characteristic is the concentration parameter, Cmr, 
defined as the ratio between the virial and inner radii, 

Cvir — Rvir/'^s- (3.13) 

A third relation between the parameters of the NFW profile is 

M^ir = ATTpsrlA{Cmr), A{Cmr) = ln(l + Cmr) " —^^H—. (3.14) 

J- -f- Cyif 

The three relations (Eqs. 13.121 13.131 and I3.14p allow the usage of any pair 
out of the parameters defined so far as the two independent parameters that 
fully characterize the profile. Finally the circular velocity curve for the halo 
is translated by 

where x = r/rg. The maximum velocity occurs at a radius Vmax — 2.16rs 
and is given by l^^ax/Kfir - 0-216 c^,^r•M(c«r)• 

Although in principle the NFW is a two-parameters family of curves, 
from statistical analysis of the simulated halos it turns out that there is an 
anti-correlation between the concentration and the halo mass |88]. Following 
[891 at z = one obtains: 
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that leads to rg ~ 27 yjQTT^j kpc- The concentration Cmr increases with 
the redshift of formation while decreasing with the halo mass, thus fulfilling 
the hierarchical clustering requirements. 

In Fig. 13.111 the NFW circular velocities from Eq. 13.151 are shown, using 
the relation 13.161 for different values of the virial mass. Each curve cor- 
responds to Mjjir = 1O^^1O"/^M0, with n = 1 . . .9 from the lowest to the 
highest curve. It is also shown the NFW maximum velocity dependence with 
the same virial mass range. 
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Figure 3.11: The NFW circular velocities with the radial coordinate in phys- 
ical units. Each curve corresponds to M^r = lO^^lO'^/^M©, with n = 1 . . . 9 
from the lowest to the highest curve. Also shown the NFW maximum veloc- 
ity dependence with virial mass. 

The NFW result has been confirmed by a number of subsequent studies 
(see e.g. [201 ED E21 [MI although there is some disagreement regarding 
the innermost value of the logarithmic slope 7. NFW argued that a fitting 
formula where 7 = {l + 3y)/{l + y) (where y = r/r^ is the radial coordinate 
in units of a suitably defined scale-radius r^) provides a very good fit to the 
density profiles of simulated halos over two decades in radius. Some authors 
(see [SHI im [96]) have argued that 7 converges to a value of ~ —1.5 near the 
center, rather than —1 as expected from the NFW fit. Others [9^ initially 
obtained much shallower inner slopes (7 ~ —0.7) in their numerical Simula- 
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tions, but have now revised their conclusions; these authors now argue that 
CDM halos have steeply divergent density profiles but, depending on evo- 
lutionary details, the slope of a galaxy-sized halo at the innermost resolved 
radius may vary between —1.0 and —1.5. 

3.6 The cusp vs core issue 

Although the existence of DM has been inferred for several decades, it is only 
recently that we start to shed light on crucial aspects of the DM distribution. 
Initially, the main focus was on the presence of a dark component [51]; this 
later shifted to investigating the ratio of dark to visible matter [53. 98l [99] . 
Today, the focus is mainly on the actual density profile of dark halos (e.g. 

Any successful cosmological model must be able to reproduce both ob- 
served large and small scale structures, from galaxy clusters to galaxy halos. 
A fundamental prediction of the cosmological CDM simulations is that viri- 
alized DM halos have an universal spherically averaged cuspy NFW density 
profile that disagrees with a number of observations. Such cusps in the DM 
distribution would certainly have very interesting implication for particle 
DM searches. For example, it could be possible to detect gamma rays from 
annihilations of very heavy DM particles in the centre of our Galaxy (e.g. 
|102] ). and present limits on radio and gamma-ray emission from the Galac- 
tic centre would then significantly constrain the mass of DM particles such 
as neutralinos. This section is devoted to address the cusp vs core issue, 
that has stimulated a lot of discussions as it has the potential to provide 
interesting new insights into the nature of DM and its possible interactions 
with visible matter (for reviews, see [1031 1104] ). 

A cored distribution, i.e. a density profile flat out to a radius that is 
a significant part of the disk size, has been often adopted (and represented 
by an isothermal profile, e.g. [105] ). although the implications of this dis- 
tribution appeared only after that cosmological A^-body simulations found 
that CDM virialized halos achieve a cuspy density profile. When the first 
simulations of CDM halos became available (e.g. [106] ). they had a central 
density profile approximately p{r) oc r~^, which has come to be known as 
the central "cusp". 

The structure of the inner regions of galactic halos was soon investigated 
by [1021 1107] . who used RCs measurements of some DM dominated dwarf 
galaxies (see Figs. 13.121 and l3.13p . It was pointed out a tension between the 
kinematical data and the predictions of simulations: DM halos seemed to 
prefer cored density distributions rather than cuspy ones. 

To cope with this observational evidences, [108] proposed an empirical 
profile (see Eq. 13.91) that successfully fitted the halo of those RCs, the so- 
called "Burkert profile" and since then has been mostly used to represent 
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Figure 3.12: DM contribution to the circular velocity of two dwarfs, as a 
function of distance from the center in units of the Hl-disk scale length 
|102| . Lines show the radial behaviour assuming a DM Hernquist (dotted), 
r~^-^ (dashed), and constant (solid) density profiles. 

cored dark halos. 

Meanwhile, theorists have done simulations with increasing resolution. 
On the basis of simulations with tens of thousands of particles per DM halo, 
NFW [U [84] showed that halos from galaxy to cluster scales have density 
profiles that are described fairly well by the fitting Eq. 13.111 

An extensively 'galaxy by galaxy' comparison then started between the 
predicted NFW density distribution and those actually detected for the dark 
halos around disk galaxies highlighting a CDM crisis and becoming the main 
goal of several pubhcations d El EH HOOl \mi \lM EH \ll2i [IH flTHl 
[TTHlfTTTt fTTH]. 

In the work of [11911120] an excellent sample of high-quality optical RCs, 
satisfying the following quality requirements were used to obtain the veloc- 
ity profiles of the surrounding DM halos: i) data extend at least out to 
the optical radius, ii) they are smooth and symmetric, iii) they have small 
internal rms, iv) they have high spatial resolution and a homogeneous ra- 
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Figure 3.13: Observed RCs as a function of galactocentric distance [107| . 
Lines show an approximately isothermal dark halo fits {p{r) oc l/(r^ +''^), 
where Vc is the core radius) to the RCs before (dotted) and after (solid) in- 
cluding the luminous contributions. Dashed lines are obtained with a Hern- 
quist profile. 



dial data coverage of 30-100 data points between the two arms, v) each RC 
has 7-15 velocity points inside Ropt, each one being the average of 2-6 inde- 
pendent data, vi) the RCs spatial resolution is better than 1/20 Ropt, the 
velocity rms is about 3% and the RCs logarithimc derivative is generally 
known within 0.05. It was found that they increase with linearly with radius 
at least out to the edge of the stellar disk, implying that, over the entire 
stellar region, the density of the dark halo is about constant. The mass dis- 
tribution was modeled as the sum of a stellar Freeman disk and a spherical 
halo (y^ = + ^J")' whose contribution to the circular velocity is given by 
[Ml M- Vni^) = - + a^)xy{x^ + a^), where x = r/Ropt, a is 

the core radius measured in units of Ropt and /3 = {V^/V'^)^^^^. It has been 
shown (e.g. [1211 1122] ) that by taking into account the logarithmic gradient 
of the circular velocity field defined as: V(r) = ^aiogr^^ ' significantly 
increase the amount of information available from kinematics and stored in 
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Figure 3.14: URC fits (thick solid line) to the RCs (points with errorbars) 
within the Constant Density Region |120| . Thin solid lines represent the disk 
and halo contributions. The maximum disk and the minimum disk solutions 
are also plotted (dashed lines). 



the shape of the RC. was calculated on both velocities and logarithmic 

rrrQrl^or,^-o■ a/2 — Vi-Vmodeliri;P,a) i 2 _ V^n-v V(n)- Vmodei ;/3,a) 

graaiems. Xv — 2^i=i Wi ^"^^ ~ Z^i=i W] ' 

and the parameters of the mass models derived by minimizing a total Xtot^ 
defined as: Xtot = Xy + Xy- The derived mass models are shown in Fig. 
13.141 alongside with the separate disk and halo contribution. It is clear that 
the halo curve is steadily increasing out to the last data point. Note also 
the uniqueness of the resulting halo velocity model: the maximum-disk and 
minimum-disk models almost coincide. 

This work is complementary to that of pUJO] who derived for 140 objects 
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Figure 3.15: The dark halo slopes Vh as a function of Vopt [100]. As a 
comparison, in CDM -0.1 <Vh < 0.5. 

of different luminosity Vh, the logarithmic gradient of the halo velocity at 
Ropt (see blue points in Fig. 13.151 red points represent the results from 
[1]). The results are impressive: the halo mass profiles at Ropt turn out to 
be i) independent of the galaxy properties, ii) Universal and in) essentially 
featureless in the sense that for any spiral the stellar disk is embedded within 
a constant density sphere. 

The highest possible value for v^^^"^ is 0.5, that is achieved on the 10 
kpc scale only for c < 5 (see [88[I123| ). i.e. for low values of the concentration 
parameter, a property of low-O universes. This value is quite inconsistent 
with the average value found in spiral dark halos, especially if one considers 
that high resolution N-body simulations converge to a maximum value of 



Of crucial importance is also the absence of a significant scatter in the 
Vh vs. logVopt relationship. In fact, the CDM theory predicts that, in a 
very wide region centered at ~ 10 kpc and including Ropt independently of 
its relation with the virial radius, galactic halos with the same mass do not 
follow a unique velocity curve but a family of them. These can be described 
by a set of straight-lines with slopes varying between —0.1 and +0.5 (e.g. see 
Fig. 6 of [BSj- According to CDM the V h — logVopt plane should be filled well 
beyond the tiny strip of Fig. I3.15[ Taken at its face value, the observational 
constraint variance [Vh < 0.1) could imply, within the CDM scenario, that 
protospiral halos are coeval and have similar merging histories. A second 
possibility may be that the disk length-scale Ropt, in units of virial radius, 
is strongly coupled with the structure of the DM halo (e.g., [1241 11251 1126] 
but see also [88]). 

Fig. 13.161 shows the URC and NFW halo fits to the RCs, leaving c and 
Vs as free parameters, constraining a conservative halo mass upper limit of 




1/4 p] 
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Figure 3.16: NFW best fits (solid lines) of the RCs (filled circles) compared 
with the CDR fits (dashed lines) |120| . The values are also indicated. 

2 X W^'^Mq: for most objects the NFW models are unacceptably worse than 
the URC solutions, and the resulting CDM stellar mass-to-light ratios turn 
out to be in some cases unacceptable low. See the particular case of the ESO 
116-G12 galaxy in Fig iXTTl 

Particular attention has been extensively given to RCs of Low Surface 
Brightness galaxies, where the implied cosmological parameters from NFW 
halos are inconsistent with ACDM picture, in what the observed concen- 
trations of the NFW halos are too low |127l 11281 1129| . Furthermore, much 
better fits to LSB observations are found when using cored halo models 

[9| [ml [TT21 M fTM [ml [T321 m [TM m 

Fig. 13.181 plots the derived mass profiles of the high-resolution LSB RCs 
sample of |112] . It is clear that most of the galaxies are characterized by an 
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Figure 3.17: Left: ESO 116-G12 density dark halo [M]. Right: CDM pre- 
diction. 

almost flat inner core with a radius of a few kpc, in contrast with the steep 
inner a = —1.0 power- law slope of the NFW profile. The values of the inner 
slope are plotted against the value of in Fig. 13.191 showing that these 
galaxies are consistent with cored halos. 

The analysis of [133| on high resolution Ha and HI RCs of 4 late-type 
dwarf galaxies and 2 LSB galaxies, based on different halo models, is shown 
is Fig. I3.20t their findings are in favour of a Burkert profile. NFW and the 
Moore profiles are inconsistent with the observed RCs in the inner regions 
in what they both predict a too fast rising RC because of the presence of the 
cuspy cores. 

Owing to the many steps in the data analysis, however, there could be 
subtle systematics errors that could distort the results, or in any case render 
the results poorly constrained. This has triggered the debate concerning 
the reliability of the data and the question of how well the mass models are 
constrained. 

The earliest observations which indicated cores in LSB galaxies were two- 
dimensional 21 cm HI velocity fields |102l llOTl I137| . Beam smearing (i.e., 
low spatial resolution) was suggested to be a systematic effect that would 
erroneously indicate cores |138l 1139] . This question was addressed by long- 
slit Ha observations which had an order of magnitude increase in spatial 
resolution (see, e.g., [1111 1140] ): yet cusps did not appear, showing that 
beam smearing had been of only minor importance in the HI observations. 
Possible systematic errors in the long-slit spectroscopy (e.g. [1181 11411 1142| 
like slit misplacement [143] and non-circular motions have since become the 
main concern. 
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Figure 3.18: Mass profiles derived from HI high-resolution LSB RCs |112| . 
The profiles can be characterized by a steep outer component, and a 
shallower inner cored component. Lines represent the best-fitting minimum 
disk models: pseudo-isothermal (full), the NFW (long-dashed) and a power- 
law fitted to the inner shallow part (thick short-dashed). The slope a is 
given in the top-left corners of the panels. The arrows indicate an angular 
size of 2", the typical value of the seeing. 



An extensive modeling was then conducted in which the RCs of both 
cuspy and cored halos were subjected to various effects and the conclusions 
were that no systematic effect will entirely mask the presence of a cuspy halo 
for realistic observing conditions (see |132| . but also |143| ). 

There are also claims that the observations could actually be consistent 
with the DM density profiles predicted by the CDM simulations, not only 
by considering the Hi data alone |138l 1144] . but also by combining Ha and 
HI data [1431 1145| . This is the reason why particular care should be taken 
in choosing a suited sample and in performing the data analysis. Note that 
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Figure 3.19: Value of the inner slope a of the LSB mass-density profiles 
plotted against the radius of the innermost point |112| . Over-plotted are 
the theoretical slopes of a pseudo-isothermal halo model (dotted lines) with 
core radii of 0.5 (left-most), 1 (centre) and 2 (right-most) kpc. The full line 
represents a NFW model, the dashed line a CDM r~^'^ model. 

recent simulations (e.g. |146| ) do not converge to a well-defined value of 
the inner slope down to the resolution limit, even though the slope of the 
DM density profile (defined as —dlnp/dlnr) at 1% of the virial radius is still 
about 1.2 for a typical galaxy. Notice also that the observational results on 
spiral galaxies show a discrepancy with the standard ACDM predictions well 
beyond the resolution limit of the simulations. 

In view of these discrepancies, together with the missing satellite problem 
(see e.g. [T2l 11471 1148j ). many alternatives to the CDM paradigm have been 
proposed. These include broken scale-invariance [1491 1150] . warm DM [1511 
I152| . scalar field DM [153111541 1155^[T56] . and various sorts of self-interacting 
or annihilating DM [HZl UMl \Mi [Ml \Mi USSl \M\- Whereas particle 
physics does not prefer CDM over these alternatives, it has the advantage 
of having no free parameters. Furthermore, most of these alternatives seem 
unable to solve both problems simultaneously [1521 11641 11651 [T66] . and face 
their own problems \T67\ [T68l [T69l [TTOl [TTT] . On the other hand, there 
are claims that the sub-structure and core problems might be solved once 
additional baryonic physics are taken into account. Several studies have 
suggested that processes such as reionization and supernova feedback can 
help to suppress star formation and to decrease central densities in low-mass 
DM halos (e.g., [T72l [T73l [T38l [T7il [T75] ). 
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Figure 3.20: RCs of the LSB sample of [I33j compared to models in minimum 
disk hypothesis: asterisks are the Ha points, open squares are the HI data; 
the continuous line is the Burkert profile, the dashed line is the NFW profile, 
and the dotted-dashed line is the Moore profile; the numbers in parenthesis 
are the Xu the three models 



Whereas these processes may indeed help to solve the problem with the 
over- abundance of satellite galaxies, the suggestion that feedback processes 
can actually destroy steep central cusps seems somewhat contrived in light 
of more detailed simulations. For instance, as shown by |172| . the effects are 
only substantial if large fractions of baryonic mass are expelled, which seems 
hard to reconcile with the low ejection efficiencies found in more detailed 
hydro-dynamical simulations (e.g., |1761I177] ). In the recent work of [178], it 
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Figure 3.21: DDO 47 RC (filled circles) best-fitted by Burkert halo + stellar 
disk (solid line) and by NFW halo + stellar disk (dashed line) mass models 



is proposed that angular momentum transfer from baryons to DM during the 
early stages of galaxy formation can actually flatten the halo inner density 
profile and modify the halo dynamics. 

It is evident from the above discussion that the long-time popular CDM 
paradigm is currently facing its biggest challenge to date. The literature on 
the cusp VS core issue is vast and there is a clear consensus on that up to 
now there is not even one spiral that requires a NFW halo. 

The importance of the issue, that concerns the very nature of DM, and 
the fact that these early results were questioned on several different aspects, 
has triggered new investigations characterized by the study of few proper 
test-cases with higher quality kinematical data, by means of properly devised 
analysis [6j. These improvements were absolutely necessary in what to obtain 
reliable DM profiles requires extended, regular, homogeneous RCs reliable up 
to their second derivative and free from deviations from the axial symmetry. 
Then, up to now, few tenths of objects have qualified to undergo such critics- 
free investigation (e.g. the list in [61171 11301 11791 ITOj ). In all these cases data 
and simulations were found in plain disagreement on different aspects: the 
best-fit disk + NFW halo mass model 

• fits the RC poorly and it implies 

• an implausibly low stellar mass-to-light ratio and 

• an unphysical high halo mass. 
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Figure 3.22: Halo pseudoisothermal (ISO) and NFW fits to the RCs sample 
of 1^ (circles), together with previous Ha (stars and squares) and HI data 
(triangles). N FWconstrained refers to fits that are constrained to matching 
the velocities at the outer radii while constraining reasonable values for the 
concentration parameter. 



As an example, it is worth to discuss in detail the case of the nearby 
dwarf galaxy DDO 47 |179| . The HI observations have adequate resolu- 
tion and sensitivity, showing that the HI 2D kinematics is very regular. 
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Figure 3.23: DM density profiles for DDO 47 and ESO 287-G13 (dots), as 
yielded by the best fits Burkert halo [180j . Solid lines: NFW density profile 
such that the mass inside the last measured point is equal for the two profiles. 
Dashed line: best-fit NFW for ESO 287-G13. 



with a well-behaved velocity field. The observed velocity along the line of 
sight V|os has been decomposed in terms of harmonic coefficients: Vjos = 
Co + '^^=i[cj cos{ jip) + Sjsm{jtp)] where ip is the azimuthal angle, cq is the 
systemic velocity, ci is the rotation velocity; it is found that the coefficients 
si) S3 32 have a small amplitude that excludes significant global elongation 
and lopsidedness of the potential and detects non-circular motions with am- 
plitude and radial profile very different from that necessary to hide a cuspy 
density distribution in the observed RC. The RC mass modeling, shown in 
Fig. 13.211 finds that the DDO 47 dark halo has a core radius of about 7 
kpc and a central density po = 1-4 x 10~^^ g cm~^, i.e. a much shallower 
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distribution than that predicted by the NFW profile. Is it possible to fit the 
data with Eq. 13.111 at all? It gives the wrong shape and the clear prediction 
of CDM is simply not realized: 

PNFW / Pobs- (3-17) 

Fig. l3.22l shows the best-fit pseudoisothermal {piso{R) = Pq/[^+{R/ Rc)'^]) 
and NFW halos to the sample of [9] of high-resolution 2D optical data com- 
bined with previous long-slit and HI RCs: the cored halo better represents 
the data in most cases. Moreover, the NFW concentrations are mostly found 
to be too low. 

An accurate mass modeling of the external regions of a couple of test-case 
spirals and a careful determination of the densities and enclosed masses of 
the DM halos at the farthest radii of 37 high quality RCs [180| . has brought 
to the discovery of a new problem for the NFW halos. In addition to the 
well-known evidence for which in the inner regions of galaxies {R < 2Rd) 
the DM halos show a flattish density profile, with amplitudes up to one order 
of magnitude lower than the ACDM predictions, at outer radii {R > 4i?D) 
the measured DM halo densities are found higher than the corresponding 
ACDM ones (see Fig. This implies an issue for ACDM that should be 

investigated in future, when, due to improved observational techniques, the 
kinematic information will be extended to the ~ 100 kpc scale [10]. This new 
discrepancy provides additional information on the nature of the cusp vs core 
issue: self-interacting or annihilating DM proposed as the cause for the inner 
discrepancy may be in difficulties in that it will cause a rapid convergence to 
the NFW profile in the luminous parts of galaxies and beyond once a critical 
density value is reached. 



3.7 Final remarks: intriguing evidences 

The distribution of luminous and DM in galaxies shows amazing properties 
and a remarkable systematics hat make it as one of the hottest cosmological 
issues. There is no doubt that this emerging observational scenario will be 
decisive in guiding how the ACDM-based theory of galaxy formation must 
evolve to meet the challenge that the observational data are posing. 

In all cases studied up to date a serious data-prediction discrepancy 
emerges, that becomes definitive when we remind that the actual ACDM 
halo profiles are steeper than the standard NFW ones considered here and 
that the baryonic adiabatic collapse has likely contracted them further. As 
a final remark I present in Fig. 13.241 a plot of the logarithmic circular ve- 
locity slope out to 6 scale lengths, of a stellar disk -|- halo model, defined as 

^ . The blue line represents the Freeman disk for a typical massive spi- 
ral; Dashed coloured (dotted) lines represent the NFW+disk model (URC), 
for a typical object with high (magenta) and low (red) luminosity. Dashed 
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black line represent NFW halo for a low luminosity object. Dot-dashed line 
represent a flat RC (V=const). Prom these simple figure it clearly comes out 
that the "flat rotation curve" paradigm is not only a completely wrong as- 
sumption, but on top of it models having this starting point loose all crucial 
information of distinguishing one model to the other. 




Figure 3.24: Logarithmic circular velocity slope as radius function in scale 
length units. Lines represent: full blue, the Freeman disk for a typical mas- 
sive spiral; dashed (dotted) coloured, the NFW-|-disk (URC) model for a 
typical object with high (magenta) and low (red) luminosity; dashed black, 
the NFW halo for a low luminosity object; dot-dashed, a flat RC (V=const). 
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Chapter 4 

Alternatives to Dark Matter 



Despite the important achievements of the CDM paradigm, as described 
in the previous chapters, this model requires however that the the overall 
dynamics of the Universe is dominated by two mysterious forms of matter 
and energy. 

As far as the only evidence of dark matter or dark energy is of gravita- 
tional origin it is reasonable to imagine that what we observe is neither a 
new form matter nor energy but a deviation of the law of gravitation from 
General Relativity. It has to be noticed that deviations from general rela- 
tivity are well motivated also from a pure theoretical point of view (string 
theory for example requires extra dimensions, possibly even "large", which 
may cause deviations from the Newtonian gravitational attraction in the 
sub-millimeter scale). 

In the this chapter two models of alternatives to dark matter are briefly 
discussed: MOND (MOdified Newtonian Dynamics) and f{TZ) theories. 

4.1 MOND 

MOND ( [181j . see |182| for a review) is certainly one of the most studied (and 
successful) model of modification of gravity. In its original formulation it 
was a pure phenomenological description (without o proper "theory" behind) 
constructed to explain two observational systematics of spiral galaxies: 

• the misleading paradigm of flat rotation curves 

• the existence of a relationship between rotational velocity and lumi- 
nosity (the Tully-Fisher relation) which implies a mass-velocity rela- 
tionship of the form M oc V", with a ~ 4. 

The basic idea of MOND is to introduce a fundamental acceleration scale, 
ao, below which deviations from the Newtonian dynamics appear; previous 
studies |183j found that ao ~ 1-2 x 10~^ cm s~^ (notice that oq ~ cHq). 
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Two formulations of the original model exist depending on whether the law 
of inertia or gravity is modified. Notice that in both cases the model clearly 
does not respect the strong equivalence principle. 

In the first case the modified law of inertia proposed by Milgrom is: 

mafi{—) = F, (4.1) 

Go 

where ^(x) is an interpolation function whose asymptotic values are fi{x) ~ x 
for X <C 1 (called Mondian regime) and fj,{x) ~ 1 for x » 1 (called Newtonian 
regime) . 

Formulated as a modification of gravity the model takes the following 
form: 

ff/i(M) = (4.2) 
ao 

where g is the effective acceleration and qn is the standard Newtonian ac- 
celeration. 

In the Mondian regime the effective gravitational acceleration takes the 
form g = ^/givao- Assuming a gravitational field generated by a point source 

2 

of mass M and imposing the condition for circular orbits {g = the circular 
velocity can be calculated: 

v'^ = GMao, (4.3) 

where G is the gravitational constant. Notice that the circular velocity does 
not depend on the radius in the Mondian regime (and hence the concept 
of fiat RCs) and a Tully-Fisher relation is obtained respecting the original 
requirements of the model. 

Even though in a general case a modified version of the Poisson equation 
should be solved, Eq. 14.21 can be shown to be a good approximation for 
axisymmetric disks |184| . The interpolation function has been given usually 
the following functional form: 

/^orig(a;) = (4.4) 

However, it is obvious that a whole family of functions are compatible 
with the required asymptotic behaviours. For instance, |185| proposed that 

f^FBix) = (4.5) 

1 + X 

could be a better choice in what it is compatible with the relativistic theory 
of MOND put forward by Bekenstein [ES]. [I87| showed that Eq. leads 
to a slightly different value of cq: oq = 1.35 x 10~® cm s~^. 

The model as described above faces a fundamental conceptual difficulty 
in analyzing composite systems, in fact, at the microscopic level the charac- 
teristic acceleration of atoms and molecule is never in the Mondian regime. 
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However although micro-systems are not Mondian the composite system in 
certain circumstances is. Recipes are then necessary for the apphcabihty of 
MOND on many body systems where the internal accelerations are above 
the MOND scale while the external acceleration is below that scale. For this 
reason MOND as conceived in its original form can not be considered as a 
satisfactory theory. 

It is clear that RCs are not asymptotically flat as originally assumed in 
construction of MOND. An asymptotically gently decreasing RC however 
can be obtained whenever the disk surface density is of order ~ oq/G or 
above. In this case in fact the internal accelerations of the disk breaks the 
Mondian regime allowing for a quasi-Newtonian decline. 

Although MOND is constructed to obtain "flat" RCs, it is able to fit 
a number of RCs and in many cases it correctly predicts general scaling 
relations linked to RCs fTMl 



[TMITMfTMfminMlfTMj . 

In the case that Eq. 14.41 is used as the interpolation function, then within 
the MOND framework the observed circular velocity lobs('') can be expressed 
as a function of oq and the Newtonian baryonic contribution Fbarl'^) to the 
RC: 

/ I / , \ 



1 + A 1 + 



_2rao_ 

V"^ (r) 

bar ^ 



1 



v 



(4.6) 



Kitars('') + ^sas(^) (ignoring the contribution of the bulge) 



where V^^^{r) 

Kitars('^) and Vgss{r) are the Newtonian contributions to the RC of the stellar 
and gaseous disks, respectively (see [181| ). The amplitude of K)tars('^) can be 
scaled according to the chosen, or fitted, stellar mass-to-light (M/L) ratio. 
^gas(^) is derived from HI observations, when they are available. 
If instead Eq. 14.51 is used the equivalent of Eq. 14.61 becomes: 




(4.7) 



(see e.g. |194| ). Note that the second term of the right-hand side of Eqs. 
14.61 and 14.71 acts as a "pseudo-dark matter halo" term and it is completely 
determined by the luminous matter. As expected, it vanishes in the limit 
ao 0. 

The MOND model can be applied also to pressure-supported systems. 
Assuming an isotropic isothermal system of total mass M, it is possible to 
obtain a relation between the radial velocity dispersion Ur and the mass 
density distribution p: 

-2 



dln{p) 
dln{r) 



(4.8) 
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Figure 4.1: On the left is a color image of the merging cluster 1E0657— 558, 
with the white bar indicating 200 kpc at the distance of the cluster |198| . 
On the right is an imagine of the same cluster from Chandra measurements. 
Shown in green contours in both panels are the weak lensing reconstruction. 
The blue +s show the location of the centers used to measure the masses of 
the plasma clouds. 

which in turn implies: 



which is similar to the observed Faber-Jackson relation for elliptical galaxies. 

Clearly the empirical formulation of MOND faces difficulties when com- 
pared with effects which originates from General Relativity (i.e. space-time 
curvature), in particular the gravitational lensing. To properly study these 
effects a covariant MOND formulation is necessary. After several attempts 
in 2004 (more than 20 years after the original MOND formulation!) it was 
proposed by J. D. Bekenstein a Tensor- Vector-Scalar (TeVeS) field theory 
which correctly reproduces the main MOND features at small accelerations 
while preserving the Newtonian regime for higher accelerations |186| . This 
theory allows not only the comparison of MOND with the lensing measure- 
ments but also the development of a theory for structure formation and a 
comparison with the modern cosmological measurements. 

Unfortunately the resulting theory is unable to fit the WMAP results in 
a pure baryonic framework and a form of DM is necessary [195| . The kind of 
DM required by MOND however may well be represented by neutrinos with 
a mass of 2eV (see e.g. [196j ). very close of the present experimental limits 



A serious challenge to the MOND-TeVeS theory came recently from the 
measurements of two merging cluster of galaxies (called "Bullet Cluster" 
|198| ). In this cluster it has been shown by means of gravitational lensing 
as well as x-rays measurements that the lensing source is misplaced with 
respect to the baryons (see Fig. I4.ip . Also in this case however the presence 
of neutrinos with a mass of 2eV would rescue MOND from a failure |196| . 




(4.9) 



m- 
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The on-going experiment KATRIN |199| however will be able soon to probe 
the range of neutrino masses necessary for MOND. 

Clearly the necessity of a form of DM also in MOND completely remove 
the beauty of the original proposal (although not excluding the model com- 
pletely). 



4.2 /(7^) theories 

This class of theories of gravitation has been studied since the very begin- 
ning of the appearance of General Relativity [2001 1201| . General Relativity 
as formulated by Einstein suffers from well known problems (among them 
the problem of singularities and the lack of a full quantization). The early 
attempts in modifying Einstein original theory went in the direction of solv- 
ing the above mentioned problems |202l 12031 1204| . More recently however 
the discovery of the dark components of the Universe modified the trend mo- 
tivating theorists in finding a gravitational theory which correctly describes 
the observed space-time geometry without the hypothesis of unknown mi- 
crophysics. 

Considering the Einstein equations: 

R,.u - ^9„uR = SttGT^^, (4.10) 

where G is the Newton Constant, R^i, and R are the Ricci Tensor and 
Scalar and T^^i, is the stress-energy tensor, the dark component is defined by 
the difference between the observed stress-energy tensor and the measured 
geometrical quantities: 

(-^Mi')darfc ~ g^Q (^Riiu ~ -^9iivR^ — iTfiu)yis^ ■ (4-11) 

The idea behind the modern version of f{R) theories is that iT^u)^ark 
nates from a modification of the fundamental Einstein equations rather than 
a new form of Dark Energy or DM. 

The Einstein equations are obtained from the Einstein-Hilbert action 
which reads: 

Seh = [ d^xV^R. (4.12) 



From this action however two variational methods can be used to obtain the 
correct equations: the metric approach (where the variation is considered on 
the metric) or the Palatini approach (where variations of the metric and of 



the connection are assumed to be independent, for reviews see [2051 1206] , 

In f{R) theories the starting action is a straightforward generalization of 
the Einstein-Hilbert one: 

^^^itcf d'^^fiR)- (4.13) 
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In this case however the two variational approaches lead to different equa- 
tions of motion. Actually for these kinds of theories a third variational 
method exists called metric affine f{R) gravity [207] where also the matter 
action (not included in Eq. I4.12p is supposed to depend on the connection. 
In what follows only the metric f{R) theory will be considered. 

In general the action will be the sum of the contribution from the theory 
of gravity and the contribution from the theory of matter: Stot = + Sm- 
Applying the metric variation the following equations are obtained: 

f'{R)Rf,u - \f{R)g^u - [V^V, - g^M f'{R) = kT^,, (4.14) 

where T^,y is the usual stress-energy tensor of the matter action. Clearly 
when f{R) = R the usual Einstein equations are obtained. The idea behind 
/(i?) theories is that the extra terms obtained in Eq. 14.141 is responsible for 
the dark component of the Universe. A remarkable feature of f{R) gravity 
is that it is equivalent to the Jordan-Brans-Dicke scalar-tensor theory [206]. 

Taking as an example f{R) oc i?" and assuming a FLRW cosmology the 
extra terms of the equations give rise to a term which can be cast in the 
form of the contribution from a perfect fluid with state equation of the form 
Pf =w pf [208] where: 

6n? — In — 1 
^^ = - 6n^-9n + 3 - ^^'^^^ 

Assuming certain values of n the Dark Energy value w ~ — 1 is obtained. 
Important limits on the value of n can be obtained both from cosmological 
analysis [2091 1210^ 12111 1212] and from the dynamics of the solar systems 
[2T31I2TI] 

The general treatment of f[R) theories is rather difficult and most often 
important results are obtained for specific forms of the / function. In general 
however the viability of a specific form of f[R) needs to fulfill the following 
criteria: 

• the correct Newtonian or Post-Newtonian limit must be obtained in 
the weak-field approximation 

• the theory must be stable at the classical and semi-classical level 

• the theory must correctly describe the dynamics of the cosmological 
perturbations. 



Chapter 5 



Tests for 
models 



dark matter mass 



As described in the previous chapter, there are several proposals for the 
solution for the cusp vs core controversy. In this Chapter I show how they 
can be validated or ruled out by a systematic comparison of their prediction 
with precision measurements of RCs. As an example I develop a test of one 
of these suggestions, the Gravitational Suppression model (GraS), that can 
be easily extended to the other proposals in the literature. 

5.1 Introduction 

The gravitational suppression hypothesis |215| is a phenomenological model 
that addresses the complex understanding of the DM distribution on small, 
subgalactic scales. High-resolution radio observations from spiral galaxies, 
along with their optical RCs, suggest that the DM is distributed in spherical 
halos with nearly constant density cores (see, e.g., [U [Sj [7] and references 
therein). On the other hand, theoretical predictions from the well-known 
N-body ACDM simulations (e.g., ^) present a steep density distribution 
profile in the centre of the halos: 



Tg is a scale radius and ps its characteristic density, in principle independent, 
but found related within a reasonable scatter through the halo mass, by 
the Bullock et al. [88| equation: c = Rmr/rg ~ ^'^{]^tm^)~^'^^ ^ where c 
is a concentration parameter and Rmr and Myir are the virial radius and 
mass. Mass models with a NFW density profile, given in Eq.(l), have two 
serious kinds of difficulty in reproducing the observed RCs: a) the fit is not 
satisfactory, i.e., x^ed ^ (^^^j ^-S-; [SI references therein); b) the values 
of the parameters of the best-fit mass models are clearly unphysical. In detail. 



Phalo{r) 



Ps 



(5.1) 



(r/rs)(l +r/rsY' 
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the values for the halo mass result much higher than those we obtain from 
weak lensing halo models [216] and from the analysis of galaxy baryonic mass 
function [77]: Mhaio ~ 3 x W^^Mq {Lb/W^^Lq)^/'^. In the same way the 
values of the disk mass-to-light ratio result much lower than those derived 
from colours of spirals [77l[2T7l[2T8| : log(MD/M0) « -1.6-M.21og(LB/L0), 
i.e., 0.7 < Md/Lb < 4. 

Several solutions have been proposed for the above issue, most of them 
related either to a better comprehension of structure formation (e.g., |219| ) 
or to new fundamental physics (e.g., [158] ). Alternatively, the presence of 
noncircular motions in galaxies has been advocated to reconcile (up tp 70% 
in the Low Surface Brightness of) the observed kinematics with the cuspy 
density profile (e.g., [2201 [M], but see also [T79] l. 

5.2 The Gravitational Suppression model 

The original proposal by Piazza & Marinoni (PM) GraS model, instead, 
modifies the usual Newtonian potential of the DM felt by baryonic test par- 
ticles in such a way that the NFW kinematics and the observed one become 
in agreement. According to PM, the NFW profile is used because GraS does 
not affect the DM dynamics, but only the dynamics in the mixed sector 
DM-baryons, so both primordial DM perturbations and halo formation are 
unaffected, and well-known N-body simulation results can be assumed. The 
idea is adding a Yukawa contribution to the gravitational potential 

^'^4>Newton = 47rG {pbaryons + Phalo), (5.2) 

from a hypothetical short-range interaction just between dark and luminous 
matter 

(V^ - A"^) 4>Yukawa = 47rG phalo-, (5-3) 

where A is a scale range parameter. The effect is damping the gravitational 
interaction on small scales. The final potential is then 

i'halo = 4>Newton + « (/^Yukawa- (5.4) 

a is a strength parameter and taken to be —1 in order to have the maximum 
possible gravitational suppression |222| . The circular velocity is related to 
the potential by 

^halo = ^halo, Newton + ^halo, Yukawa = \d4>halo/ dr\. (5.5) 

In PM model, for a (small) sample of RCs of Low Surface Brightness galaxies 
GraS was able to eliminate the above core versus cusp discrepancy. However, 
in order to allow a simple analytic calculation, they have taken a number 
of assumptions and approximations. In detail, the contribution to the grav- 
itational potential from baryons (stars and HI disk) was neglected and the 
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DM distribution was modeled with the simple form Phaio{x) = pox~^, rather 
than by Eq. (1). Further support to GraS was given in |223| where the dis- 
persion velocity of two spheroidal dwarfs (Fornax and Draco) were studied 
in this scenario. However, both large errors in the kinematic measurements 
and large geometric and orbital uncertainties of the employed mass model, 
limited the relevance of their findings. 

5.3 Data and methodology of the test 

In the present analysis of GraS we abandon the above approximations and 
test a wider and fairer sample of spirals. An in-depth review of the GraS 
model is beyond the scope of this work. Our goal is to perform a check of 
GraS. First, we assume the exact NFW profile. Second, we consider the 
baryonic contribution, so that the total potential is 

(pmodel = <t^halo + 4>disk + (j^gas, (5-6) 

where the sum of the last two terms is (pbaryons- This leads to 

^ model ^ halo ' ^ disk ' ^ gas' \'^-') 

Finally, we use a sample of high-resolution RCs of Low and High Surface 
Brightness galaxies, in order to investigate the consistency and universality 
of the model. 

Our sample represents the best available RCs to study the mass distri- 
bution of DM and it has been used in works concerning the core versus cusp 
discrepancy controversy [SI I224j . The sample includes nearby Low and High 
Surface Brightness galaxies, all poorly fitted by mass models with NFW ha- 
los that also have unphysical values for their best-fit mass parameters: DDO 
47 [US]; ESO 116-G12, ESO 79-G14 [6]; NGC 6822 [U^; UGC 8017, UGC 
10981, UGC 11455 [225]; M 31 [224]. Let us notice that in some cases H^ and 
HI RCs are both available and they agree well where they coexist. Moreover 
the RCs we analyse are smooth, symmetric and extended to large radii. 

We decompose the total circular velocity into stellar, gaseous and halo 
contributions, according to Eqs.(l)-(5), where the latter contains the addi- 
tional DM-baryons interaction. Available photometry shows that the stars 
in our sample of galaxies are distributed in a thin disk, with exponential 
surface density profile S£)(r) = (M£)/27ri?|,) e~~^l^° ^ where Md is the disk 
mass and Rr, is the scale length. The circular velocity contribution is given 
by V;i^fc(r) = {GMd/2Rd) x^B{x/2), where x = r/Ro and G is the gravi- 
tational constant. The quantity B = IqKq — Ii/i is a combination of Bessel 
functions [63]. The contribution of the gaseous disk is directly derived from 
the HI surface density distribution. 

In a first step, the RCs are best-fitted with the following free parame- 
ters: disk mass, NFW scale radius and characteristic density, and scale range 
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of GraS. Then we redo the analysis fixing the GraS scale range parameter at 
the mean value found of A = 3.1 kpc. Notice that the published mean value 
of PM for A is quite different from ours as an effect of their simplifications: 
A = 1.1 kpc. Our value is the most favourable for the PM model: different 
values of A leads to worse performance. 

5.4 Results 

The test goes against the GraS model. For the RCs of our sample the 
NFW mass halo model fails to reproduce data according to the usual pattern 
explained in the introduction. Data, not surprisingly, points to DM halos 
having inner density cores. Applying a Yukawa potential to the cuspy NFW 
halo does not solve this discrepancy. The cusp is erased and RCs are fitted 
very well, but this success is illusory in that the corresponding values of the 
parameters of the best-fit mass model remain unphysical. In table I we show 
the results of the test. We give: the values of the parameters of the mass 
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Figure 5.1: Galaxies in which GraS eliminates the core versus cusp discrep- 
ancy controversy. Y axis is the velocity in km/s. The solid line represents the 
best-fit mass model, the long-dashed line is the contribution of the DM halo, 
and the dotted and short-dashed lines are those of the stellar and gaseous 
disks. Below the RCs, we plot the residuals {Vobs — Vmodei)- 
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model and global properties of the galaxies, xled calculated with average 
typical velocity errors. In bold, unphysical values for halo mass and mass- 
to-light ratio, and xled > 2. 5. The critical density of the Universe today is 
taken to be Pcrit,o = W~^^g/cm'^. 

In detail, in the cases of ESO 116-G12 and UGC 10891, we have that 
GraS fits sufficiently well the RCs unlike the NFW, confirming that this 
model could work in some objects (see Fig. 1, table I). 

However, in the other cases, although the fits are satisfactory, the best- 
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Figure 5.2: Galaxies in which GraS does not solve the core versus cusp 
discrepancy controversy. The fitting values of the mass-to-light ratio (NGC 
6822, ESO 79-G14, UGC 11455) and halo mass (DDO 47, UGC 8017) result 
unphysical. See Fig. 1 and table I for details. 
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Figure 5.3: NFW+GraS mass model. Left: with typical values for the halo 
mass and the mass-to-light ratio. Right: with the Bullock et al. relation. 
See Figs. 1 and 2 and table I for details and comparison. 

fit values of the halo mass and mass-to-light ratio are unphysical. In fact, 
we expect (see above) the mass-to-light ratios for NGC 6822, ESO 79-G14, 
UGC 11455, to be equal to (1, 2.6, 3.5), while we found much smaller best-fit 
values (<0.02, 0.3, <0.2). In the same way, we expect halo masses for DDO 
47 and UGC 8017 to be equal to {^xIQ^^Mq, 1.9x IQI^M©), while we found 
much bigger best-fit values (8.1 x 10^^ M©, 1.5 x 10^^ Mq). Furthermore, in 
M 31 the GraS modification is negligible and irrelevant (see Fig. 2, table I). 

Let us notice that by constraining the values for the mass parameters 
within physically acceptable values, we obtain unacceptable fits for the GraS 
mass model, similar to those of the Newtonian NFW case. As an example, 
in UGC 8017 with Mhaio = 3 x lO^^M© and Md/Lb = 3Mq/Lq, GraS 
shows an unacceptable fit to data (see Fig. 3). More in general, we realize 
that for all six objects, all values of ps and within their la uncertainties 
imply unphysical halo masses and/or mass-to-light ratios. 

We now implement the Bullock et al. concentration vs halo mass relation, 
that eliminates one parameter in the original NFW profile. With this relation 
built in, GraS performs even worse than before. See in Fig. 3 the case for 
DDO 47. 

5.5 Conclusions 

In conclusion, the GraS-PM model fails to rescue the NFW profiles in a 
number of high quality well-suited RCs. Moreover, let us point out that 
there is not a pattern of this inability, so that it is presently difficult to 
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understand how to modify it in order to reach its original goal. Then the 
GraS model is a rather strong hypothesis that does not seem solve the core 
versus cusp discrepancy problem of the mass distribution of the center of 
DM halos. 

Finally, let us remark that also in this work it has emerged that the avail- 
able kinematics of galaxies is very constraining for non-Newtonian theories 
of gravity. 
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Table 5.1: Parameters of the mass models. 
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Chapter 6 



Analysis of Rotation Curves in 
the framework of gravity 

As seen in Chapter 4, modifications of the law of gravity are appeahng 
alternatives to the yet undetected DM particles. In this chapter I develop 
a test of i?" gravity in galactic scales that represents a step forward on the 
issue in what for the first time a complete analysis of a devised sample of 
Rotation Curves has been performed. 

6.1 Introduction 

It is well-known that the RCs of spiral galaxies show a non-Keplerian circu- 
lar velocity profile which cannot be explained by considering a Newtonian 
gravitational potential generated by the baryonic matter ^62j. Current pos- 
sible explanation of this controversy includes, among others, the postulate 
of a new yet not detected state of matter, the DM [226| . a phenomenological 
modification of the Newtonian dynamics [1811 12271 11821 [228] . and higher or- 
der gravitational theories (originally devoted to solve the dark energy issue, 
see e.g., [2291 [230]). 

The recent theory proposed by Capozziello, Cardone & Troisi 2007 (here- 
after CCT, |231] ). modifies the usual Newtonian potential generated by bary- 
onic matter in such a way that the predicted galaxy kinematics and the ob- 
served one have a much better agreement. They consider power-law fourth 
ordei0 theories of gravity obtained by replacing in the gravity action the 
Ricci scalar R with a function f{R) oc i?"^, where n is a slope parameter. 
The idea is that the Newtonian potential generated by a point-like source 
gets modified in to 

Cm 1 

m = - — {l + ^[{r/r,f-l]}, (6.1) 

^The term comes from the fact that the generalized Einstein equations contain fourth 
order derivatives of the metric. 
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where /3 is a function of the slope n, and rc is a scale length parameter. 
It turns out that in this theory /3 is a universal constant while Tc depends 
on the particular gravitating system being studied. In a virialized system 
the circular velocity is related to the derivative of the potential through 
V^ = r d(j){r)/dr. It is clear that (l6ll|) may help in the explanation of the 
circular velocity observed in spirals. 

We remark that any proposed solution to the galaxy RC phenomenon 
must not only fit well the kinematics but, equally important, also have best- 
fit values of the mass model parameters that are consistent with well studied 
global properties of galaxies. 

For a sample of 15 Low Surface Brightness galaxies the model described 
in CCT was fairly able to fit the RCs. However, in our view, the relevance 
of their finding is limited by the following considerations: 

• the sample contains several objects whose RCs are not smooth, sym- 
metric and extended to large radii 

• the sample contains only Low Surface Brightness galaxies while a wider 
sample is desirable 

• the universal parameter n is not estimated by the analysis itself but it 
is taken from other observations. 

In the present work we generalize the results of CCT and test a wider 
and fairer sample of spirals, improving the analysis methodology. Our goal 
is to perform a check of their model on galactic scales in order to investigate 
its consistency and universality. 

The plan of this work is the following: in Sect. 2 we briefiy summarize 
the main theoretical results described in CCT relevant for the analysis of 
our sample. In Sect. 3 we present our sample and methodology of analysis. 
In Sect. 4 the results are presented and finally the conclusions in Sect. 5. 

6.2 Newtonian limit of f{R) gravity 

The theory proposed by CCT is an example of f{R) theory of gravity [2321 
I233| . In these theories the gravitational action is defined to be: 



where g is the metric determinant, R is the Ricci scalar and Cm is the matter 
Lagrangian. They consider: 



where /o is a constant to give correct dimensions to the action and n is the 
slope parameter. The modified Einstein equation is obtained by varying the 
action with respect to the metric components. 




(6.2) 



f{R) = kR 



(6.3) 
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Solving the vacuum field equations for a Schwarzschild-like metric in the 
Newtonian limit of weak gravitational fields and low velocities, the modi- 
fied gravitational potential for the case of a point-like source of mass m, is 
given by ()6.ip . where the relation between the slope parameter n and /3 (see 
detailed calculation in CCT) is given by: 



Note that for n = 1 the usual Newtonian potential is recovered. The large 
and small scale behavior of the total potential constrain the parameter /? to 
be < /3 < 1. 

The solution (16. 1|) can be generalized to extended systems with a given 
density distribution p{r) by simply writing: 



where (pNi^) represents the usual Newtonian potential and (pdr) the addi- 
tional correction. In this way, the Newtonian potential can be recuperated 
when [3 = 0. The solution for the specific density distribution relevant for 
spiral galaxies is described in the following paragraph. 

6.3 Data and Methodology of the test 

We selected two samples of galaxies: a first with 15 galaxies, called Sample 
A, that represents the best available RCs to study the mass distribution of 
luminous and/or DM, and it has been used in works concerning modifications 
of gravity and the cusp vs core controversy ^ 12241 1234| . 

This sample includes nearby galaxies of different Surface Brightness: 
DDO 47 [179]; ESO 116-G12, ESQ 287-G13, NGC 7339, NGC 1090 [6]; 
UGC 8017, UGC 10981, UGC 11455 [225]; M 31, M 33 [224]; IC 2574 [235], 
NGC 5585 [23g, NGC 6503 [23?], NGC 2403 |238j, NGC 55 [239]. This 
sample is the most suitable for a fair test of theories like the one of CCT: 

• The RCs are smooth, symmetric and extended to large radii. 

• The galaxies present a very small bulge so that it can be neglected in 
the mass model to a good approximation. 

• The luminosity profile is well measured and presents a smooth behavior 

• The data are uniform in quality up to the maximal radii of each galaxy. 



12n2 - 7n - 1 - V36n4 _^ i2n3 - 83^2 + 50n + 1 
6n2 - 4n + 2 



(6.4) 




(6.5) 
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Let us notice that in some of these galaxies Hq, and HI RCs are both available 
and in these cases they agree well where they coexist. 

We also considered a second sample called Sample B consisting of 15 
selected objects from Sanders & McGaugh 2002 that has been used to test 
MOND. This sample consists of the following galaxies: UGC 6399, UGC 
6983, UGC 6917, NGC 3972, NGC 4085, NGC 4183, NGC 3917, NGC 3949, 
NGC 4217, NGC 3877, NGC 4157, NGC 3953, NGC 4100 [2401 [241]; NGC 
300 [2l2]; UGC 128 [2ig. Although these galaxies do not fulfill ah the 
requirements of Sample A we have analyzed them for completeness sake. The 
properties of the galaxies of the two samples are listed in table 1. Notice 
that the theory of CCT requires an analysis with a sample of high quality 
galaxies, as described above, where each luminous profile plays an important 
role, whereas this is not the case in MOND. 

We decompose the total circular velocity into stellar and gaseous contri- 
butions. Available photometry and radio observations show that the stars 
and the gas in our sample of galaxies are distributed in an infinitesimal 
thin and circular symmetric disk. While the HI surface luminosity density 
distribution ^^^^(r) gives a direct measurement of the gas mass, optical ob- 
servations show that the stars have an exponential distribution: 

EDir) = iMD/27rRl)e-^/''^, (6.6) 

where M^, is the disk mass and Rd is the scale length, the latter being 
measured directly from the optical observations, while is kept as a free 
parameter of our analysis. 

The distribution of the luminous matter in spiral galaxies has to a good 
extend cylindrical symmetry, hence using cylindrical coordinates, the poten- 
tial (|6.5p reads 

<^(r) = -G / dr' r'W) / T^A^ + o "l]}- (6-7) 



E(r') is the surface density distribution of the stars, given by (16.6(1 , or 
of the gas, given by an interpolation of the HI data points up to the last 
measured point. 13 and Tc are free parameters of the theory, with the latter 
galaxy dependent. We neglected the gas contribution to the mass density for 
radii larger than the last measured point, however we checked the goodness 
of this approximation by extending the distribution with a different kind 
of decreasing smooth curves and realized that error made in the truncated 
approximation is small enough to be neglected. 

Defining k"^ = j^^^y^, we can express the distance between two points in 

cylindrical coordinates as |r — r'| = {r + {1 — k'^ cos"^ {6 / 2)) . The derivation 
of the circular velocity due to the marked term of equation (|6.7p . that we 
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call 0/3 (r), is now direct: 

r ^ cj)p{r) = tt a {(3 - 1) G I{r), (6.8) 

where the integral is defined as 

T(r) = / dr'r''- k^-l^ i:{r') F{r), (6.9) 

with .F(r) written in terms of confluent hyper-geometric function: J-{r) = 
2(r + r') 2Fi[l, 1^ + [(^' " 2)r-' + k^r] 2Fi[|, 2, k\ 

The total circular velocity is the sum of each squared contribution: 

^CCrir) = "^N^stars + ^N,gas + ^C.stars + ^C,sas (6.10) 

where the stars and gas subscripts refer to the different contributions of 
luminous matter to the total potential (|6.5|1 . The N and C subscripts refer 
to the Newtonian and the additional correction potentials. 
Let us recall that we can write 

Vstars'^{r) = {GMd/2Rd) x^B{x/2), (6.11) 

where x = rjRi)^ G is the gravitational constant and the quantity B = 
IqKq — IiKi is a combination of Bessel functions f63|. 

Galaxies UGC 8017, M 31, UGC 11455 and UGC 10981 presents a very 
small amount of gas and for this reason it has been neglected in the analysis. 
Notice that the correction to the Newtonian potential in equation (|6.ip may 
be negative and this would lead to a negative value of Vj. In Figs. 1 and 
2 however the velocities Vc are shown only in the ranges of r where their 
square are positive. 

In a first step, the RCs are best-fitted with the following free param- 
eters: the slope (/3) and the scale length (rc) of the theory, and the gas mass 
fraction {fgas) related to the disk mass simply by Md = Mgas{l- fgas)/ fgas- 
The errors for the best fit values of the free parameters are calculated at one 
standard deviation with the Xred + the results of these fits we 

get a mean value of /3 = 0.7 ± 0.25 (n ~ 2.2). In the second step we redo 
the best-fit fixing the slope parameter at /3 = 0.7 keeping as free parameters 
only rc and fgas- Notice that in a previous paper [244] . a mean value of 
f3 = 0.58 lb 0.15 (n ~ 1.7) has been obtained, perfectly compatible with our 
result. This parameter however, is well constrained from SNela observations 
to be /? = 0.87 (n ~ 3.5), also compatible with our measurements. In our 
analysis the value /? = 0.7 is the most favorable for explaining the RCs: 
different values of f3 from the one we adopt here lead to worse performance. 
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Figure 6.1: Sample A: The solid line represents the best-fit total circular 
velocity VccT- The dashed and dotted lines are the Newtonian contributions 
from the gas and the stars, while the dot-dashed represents their sum. The 
long-dashed line is the non-Newtonian contribution of the gas and the stars 
to the model. Below the RCs, we plot the residuals (Vobs — Vcct)- See table 
1 for details. 
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Figure 6.2: Sample B: Best-fit curves superimposed to the data from selected 
objects from Sanders & McGaugh 2002. See Fig. 1 for details. 

6.4 Results 

We summarize the results of our analysis in Figs. 1 and 2 and table lH. In 
general we find for all galaxies: 

^Numerical codes and data used to obtain these results can be found at the address 
http://people.sissa.it/~martins/home.html 
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• the velocity model VccT well fitting the RCs 

• acceptable values for the stellar mass-to-light ratio 

• too vast range for values of the gas fraction (0% < fg < 100%) 

• not clear comprehension for the big variation of values for the scale 
length parameter (0.005 kpc< Tc <1.53 kpc). 

The residuals of the measurements with respect to the best-fit mass model 
are in most of the cases compatible with the error-bars, see Figs. 1 and 2, 
though three galaxies show significant deviations: NGC 6503, NGC 2403 
and M 33. 

We also find acceptable values for the B-band mass-to-light ratio param- 
eter for most of the galaxies, for which we should have approximately 0.5 < 
< 6 and a positive correlation between B-luminosity (T^ = Md/Lb; 



Table 6.1: Properties and parameters of the mass model of the analyzed 
Samples (/? = 0.7). Prom left to right, the cohimns read: name of the 
galaxy, Hubble type as reported in the NED database, adopted distance in 
Mpc, B-band luminosity in 10^ Lbq, disli scale length in kpc, gas mass in 
IO^Mq until last measured point, gas fraction in %, disk mass in 10^ M©, 
scale length CCT parameter in 10—2kpc, mass-to-light ratio in Tq, and Xred- 
The galaxies are ordered from top to bottom with increasing luminosity. 



dalaxij 


Type 




Lu 


Ro 


.v.,„. 






7'.,„. 








xU 












Sample 


A 












DDO 47 


IB 


4 


0.1 


0.6 


2.2 






96±1 


0.01 


0.6 


0.1 


0.6 


IC 2674 


SABm 


3 


0.8 


1.78 


0.6 






79±12 


0.14 


1.7±0.3 


0.2 


0.8 


NGC 6585 


SABc 


6.2 


1.3 


1.26 


1.6 






68±3 


1 


3.8±0.4 


0.7 


1.4 


NGC 55 


SBm 


1.6 


4 


1.6 


1.3 






84±7 


0.24 


2.4±0.4 


0.06 


0.1 


BSO 116-G12 


SBcd 


15.3 


4.6 


1.7 


21 






60 


2.1 


6±1 


0.6 


1.2 


NGC 6603 


Sc 


6 


6 


1.74 


2.3 






18±0.7 


10.6 


21±1.4 


2.1 


18 


M 33 


Sc 


0.84 


6.7 


1.4 


3.7 






63±2 


3.3 


7.6±0.4 


0.6 


26 


NGC 7339 


SABb 


17.8 


7.3 


1.6 


6.2 






2.8±0.2 


22 


41±7 


3 


2.3 


NGC 2403 


Sc 


3.26 


8 


2.08 


4.6 






27±0.9 


12.1 


21±1.6 


1.6 


19 


M 31 


Sb 


0.78 


20 


4.6 










180±70 


163±19 


9 


3.4 


ESO 287-G13 


Sbc 


35.6 


30 


3.3 


14 






25±1 


41 


48±5 


1.4 


3.2 


NGC 1090 


Sbc 


36.4 


38 


3.4 


100 






18±1 


47 


59±4 


1.2 


0.9 


UGC 8017 


Sab 


102.7 


40 


2.1 










9.1±0.3 


1±1 


0.2 


6.2 


UGC 11466 


Sc 


76.4 


46 


6.3 










74±3 


14±1 


1.6 


6 


UGC 10981 


Sbc 


166 


120 


6.4 










460±200 


~ IQll 


3.8 


4.9 












Sample 




3 


UGC 6399 


Sm 


18.6 


1.6 


2.4 


1 






23±3 


3.3 


10±3 


2 


0.1 


NGC 300 


Scd 


1.9 


2.3 


1.7 


1.3 






39±4 


2 


5.2±1 


0.9 


0.4 


UGC 6983 


SBcd 


18.6 


4.2 


2.7 


4.1 






24±2 


13 


4.6±10 


3.1 


0.9 


UGC 6917 


SBd 


18.6 


4.4 


2.9 


2.6 






14±1 


16 


71±17 


3.6 


0.5 


UGC 128 


Sd 


60 


5.2 


6.4 


10.7 






32±5 


23 


39±11 


4.4 


0.1 


NGC 3972 


Sbc 


18.6 


6.7 


2 


1.5 






39±3 


2.5 


2.5±0.4 


0.4 


0.1 


NGC 4086 


Sc 


18.6 


6.9 


1.6 


1.3 






44±4 


1.7 


1.4±0.3 


0.3 


1 


NGC 4183 


Scd 


18.6 


9.6 


1.4 


4.9 






60±6 


3.2 


9±2.3 


0.3 


0.3 


NGC 3917 


Scd 


18.6 


11 


3.1 


2.6 






22±1.5 


9.2±0.9 


9.8±1.4 


0.8 


1 


NGC 3949 


Sbc 


18.6 


19 


1.7 


4.1 






19±2.2 


17 


22±6 


0.9 


0.3 


NGC 4217 


Sb 


18.6 


21 


2.9 


3.3 






6.1±0.7 


52 


55±16 


2.5 


0.4 


NGC 4100 


Sbc 


18.6 


26 


2.6 


4.4 






13±1.6 


28 


20±3 


1.1 


1.5 


NGC 3877 


Sc 


18.6 


27 


2.8 


1.9 






7.3±0.8 


24 


20±4 


0.9 


0.8 


NGC 4167 


Sb 


18.6 


30 


2.6 


12 






26±2.6 


33 


26±4 


1.1 


0.5 


NGC 3953 


SBbc 


18.6 


41 


3.8 


4 






2.8±0.18 


140 


190±50 


3.4 


0.8 



6.5. CONCLUSIONS 



85 



Mf) is the disk mass and -L^ is the B-band galaxy luminosity) and [218| : 
Mn{LB) ^ 3.7 X w'^ x [{^)^-^^ giLn) + 0.095{^f^^]MQ, (6.12) 

where Lio = W^^Lbq and c/(Lb) = exp[-0.87 x (/o^^ - 0.64)^]. In detail 
we find discrepancies for NGC 55, UGC 8017, NGC 3972, NGC 4085 and 
NGC 4183. Values for the scale length parameter (rc) are in general smaller 
for less massive galaxies and bigger for more massive ones. We obtained a 
Newtonian fit for UGC 10981, as shown by the exceedingly large value for 
Tc, see Fig. 1. 

The model analyzed in this work yields better results on galactic scales 
than CDM models, where in the latter these galaxies have serious problems 
like marginal fits and unreasonable values for the stellar mass-to-light ratio, 
see e.g., p[234]. 

6.5 Conclusions 

We have investigated the possibility of fitting the RCs of spirals with a 
power-low fourth order theory of gravity of CCT, without the need of DM. 
We remark the relevance of our sample that contains objects in a large range 
of luminosity and with very accurate and proper kinematic. We find in 
general a reasonable agreement, with some discrepancies, between the RCs 
and the CCT circular velocity model, encouraging further investigations from 
the theoretical point of view. 



CHAPTER 6. ROTATION CURVES k GRAVITY 



Chapter 7 



Universal scaling relations in 
the luminous and dark mass 
distributions of spirals and 
dwarfs spheroidals 



I now turn the discussion to a more in depth analysis of DM halos in galaxies 
on a wide range of galaxy luminosity. Kinematic surveys of the dSph satel- 
lites of the Millcy Way are revealing tantahzing hints about the structure 
of DM halos at the low end of the galaxy luminosity function. In brighter 
galaxies, observations and modeling of spiral galaxies suggest that their dark 
halo parameters follow a number of scaling relations. In this work, we in- 
vestigate whether the extrapolation of these relations to the dSph regime is 
consistent with the observed internal kinematics of dSphs. The negligible 
fraction (~ 10~^ — 10~^) of baryonic matter inside the optical regions of 
dSphs is consistent with the declining trend of baryon fraction with baryonic 
(and DM) mass seen in spirals. The dSph data do not currently discrimi- 
nate between cored and cusped halos, due to our lack of knowledge about 
the anisotropy of the stellar velocity distribution and the limited spatial ex- 
tent of the stellar tracers relative to the DM. Nevertheless, although the DM 
densities in dSphs are typically almost two orders of magnitude higher than 
those found in (larger) disk systems, we find that the dSph kinematics are 
consistent with their occupancy of (cored) Burkert DM halos whose core 
radii and central densities lie on the extrapolation of the scaling laws seen in 
spiral galaxies. We discuss the potential implications of this scaling relation, 
if confirmed by future observations, for understanding the nature of DM. 
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7.1 Introduction 

DM provides the gravitational potential wells in which galaxies form and 
evolve. Over the past decades, observations have provided detailed infor- 
mation about the distribution of DM within those regions of spiral galaxies 
where the baryons reside ([75] and references therein, [551 Ell EB]). Similar 
information is now also becoming available for Low Surface Brightness galax- 
ies [HI [9]. In these disk systems, the ordered rotational motions and known 
geometry of the tracers has facilitated this achievement and an intriguing 
phenomenological picture has emerged. Spiral galaxies are composed of a 
disk surrounded by a dark halo. Inside the optical regions {R < Ropt), the 
disk is almost self-gravitating in the most luminous objects but contributes 
a negligible amount to the gravitational potential at the lower end of the 
luminosity function. Mass modeling of both individual and co-added RCs 
shows: (1) cored DM halos generally provide a better fit to the observed 
data than cusped halos; more specifically, the Burkert density profile 

P^'^ = ( ^ 2, ' (7-1) 

that contains two free parameters, the core radius tq and the central halo 
density po, reproduces the available kinematical data [4l [6l [TOl [62l 1111111121 
11331 1179] . (2) When the data for spiral galaxies are modeled assuming a 
Burkert distribution for the DM and a Freeman disk for the luminous matter, 
the parameters (DM central densities, core radii, disk masses and length 
scales) are all related by a series of scaling laws [621 [751 [77] . 

In contrast to the results gathered for disk-dominated systems, our knowl- 
edge of the mass distribution in pressure-supported systems like elliptical 
galaxies is still limited (see |245| for a recent summary of the state of art). 
However, on-going observations of Local Group dwarf spheroidal galaxies 
(dSph), which occupy the faint end of the luminosity function of pressure- 
supported systems, are currently yielding crucial information about the prop- 
erties of the dark and luminous components in these objects and, in turn, 
on the underlying physical properties of DM halos (e.g. |246l 12471 1248] ). A 
number of important questions remain unanswered. These include: 

• Is the distribution of DM on galactic (i.e. kpc) scales universal? 

• Why do the dark and luminous mass distributions appear to be re- 
lated, even though baryons dominate, at most, only the inner regions 
of galaxies? 

Work on spiral galaxies performed over the past 20 years has suggested 
some answers to these questions for this Hubble type. In contrast, observa- 
tions of the internal kinematics of dSphs have only recently begun to provide 
hints of the distribution of DM in these low-luminosity systems. The dSphs 
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Figure 7.1: Comparison of the distribution of characteristic baryonic scale 
Rd versus stellar mass Ms for dSphs (points) with the corresponding relation 
in Spirals (see [62]). For the dSphs, the stellar mass is estimated from the 
V band galaxy luminosity, assuming a stellar mass-to-light ratio of unity (in 
solar units). 



are indispensable for building up an observational picture of the process of 
galaxy formation as they extend the exploration of the dark and luminous 
mass distribution in galaxies over a much wider range of Hubble type and 
luminosity. An indication of this is given in Fig. 17.11 where we show the 
relationship between a characteristic baryonic length scale (see below for 
definitions) and the stellar mass in spirals and dSphs. The figure illustrates 
the very different ranges of baryonic mass and size scale in these two classes 
of stellar system. 

The dSphs are typically at least two orders of magnitude less luminous 
than the faintest spirals, and show evidence of being DM dominated at all 
radii. They are a primary laboratory for the bottom-up theory of galaxy for- 
mation. Moreover, being predominantly old, pressure-supported, spheroidal 
systems, their evolutionary histories are significantly different from those of 
spirals, especially in the baryonic components. There is some evidence of 
universality in the global properties of the mass distribution of dSphs. |249| 
found that the variation of the mass to light ratios of dSphs with total lu- 
minosity was consistent with the hypothesis that all dSphs contain similar 
masses of dark matter interior to their stellar distributions that implies a 
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larger proportion of DM in the less luminous objects, a main global charac- 
teristic of the spiral mass distribution, well-known since |53| . 

More recent analysis |246l I250| . based on extended velocity dispersion 
profiles rather than central velocity dispersions, have generally supported 
this conclusion. In dSphs, we note that due to the limited spatial extent of 
the stellar distributions, the radial limit of kinematic observations may be 
only a small fraction of the actual DM halo size. 

Obviously, the existence of common features, or scaling laws, relating 
the structural parameters of the mass distributions of dSphs with those of 
very different stellar systems (e.g. spirals) would be of potentially great 
significance and indicative of a Grand Picture drawn by the fundamental 
physical processes in the formation and evolution of galaxies. A number of 
recent papers have studied various scaling relations between the properties 
of hot stellar systems [2511 12521 [253] ■ A common conclusion is that the 
dSphs are outliers from other spheroidal systems in terms of many of their 
properties (mass-to-light ratios, sizes, etc.). In this paper, therefore, we 
examine whether the properties of the dSphs are consistent with another 
class of stellar system, namely luminous spiral galaxies. Additionally, while 
previous works have studied the properties of the stellar distributions, or 
global mass-to-light ratios, in this paper we make a tentative first attempt 
to compare the DM halo parameters of different systems. In particular, we 
will extrapolate the picture emerging in spirals to the region of parameter 
space occupied by the dSphs, thus comparing systems across a broad span 
of galaxy global properties and morphologies. 

The outline of the paper is as follows. In Section [721 we summarize the 
observational data used in our study and describe in detail the analysis of 
the dSph data. Section 17.31 compares the properties of the dark halos of 
spiral and dSph galaxies, while Section 17.41 discusses the relations between 
the baryonic and DM properties of these systems. Section 17.51 summarizes 
our findings and speculates on the implications for the nature of DM. 

7.2 Data 

7.2.1 Spiral Galaxies 

Our aim in this paper is to test the consistency of the dSph data with 
the scaling relations seen in spiral galaxies. When the mass distribution 
in spirals is modeled using a Burkert DM halo (with parameters po and tq) 
and a Freeman stellar disk, a tight relation between pQ and tq emerges |[62] . 
Noticeably , as it can be seen in Fig. 6, we find similar po vs tq relationships 
independently of whether the mass profiles are obtained from kinematics (i. e. 
from RCs) or from gravitational lensing data or from the analysis individual 
or coadded objects. 

More in detail here, we make use and show in Fig. 6 of the values of these 
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parameters obtained in galaxies in which the profiles have been determined 
via one of (1) the Universal Rotation Curve (see Fig. 17.61 and [77]): (2) the 
analysis of weak lensing signals around spirals (see Fig. 17. 6|) : (3) the mass 
modeling of individual RCs [Bl 171 [TT]. 



7.2.2 dSph galaxies 

The study of the internal kinematics of the Milky Way dSphs has been 
revolutionized by the availability of multi-object spectrographs on 4m and 
8m-class telescopes. Large data sets comprising several hundred individual 
stellar velocities per galaxy have now been acquired for all the luminous 
dSphs surrounding the Milky Way [Ml ESQI |2Ml ESS [256l ESZl |258l 1259] . 
The volume of the currently available data is sufficient to place reliable con- 
straints on the dynamical masses interior to the stellar distributions of the 
dSphs. However, the mass profiles are less well-determined, and the velocity 
dispersion profiles alone cannot distinguish between cored and cusped ha- 
los due to the degeneracy between mass and velocity anisotropy (see, e.g. 
[25011259] ). However, [246 ] recently showed that the kinematic data in six of 
the well-studied dSphs are consistent with their occupying cored DM halos, 
under the assumptions of spherical symmetry and velocity isotropy. Further, 
[246] note that two dSphs exhibit additional features which suggest that their 
halos are not cusped. 

Before comparing the properties of dSphs with those of spiral galaxies, 
we first re-visit the DM density profiles derived in [246] for six Milky Way 
dSphs. In particular, we investigate whether the Burkert DM profile which, 
let us recall, generally reproduces the RCs data for spiral galaxies ([4], see 
also [6]), is also consistent with the observed data for dSphs. We note that 
for spiral galaxies, it has been shown that this choice of halo model is not 
prejudicial. In the region probed by the data, for appropriate values of the 
halo parameters, actually very different from those that we actually find, the 
(cored) Burkert profile could have mimicked, to a very good approximation, 
a (cusped) NFW [3|. To proceed, we would ideally require estimates of the 
Burkert parameters po and rg, as well as their associated errors, for our 
six dSphs. However, as we discuss below, an unambiguous determination of 
whether dSph halos are cored and, if so, the sizes of their core radii, is beyond 
the scope of the present paper. Instead we will investigate the consistency of 
Burkert halos with the velocity dispersion profiles of the dSphs, as published 
in [2l6]- 

For each object, we generate 1000 random realizations of the surface 
brightness profile and velocity dispersion profile by drawing values within 
the observed error bars. We fit each surface brightness profile with a Plum- 
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Name 








PO 






M(fls3/2) 




(kpc) 


(km/s) 


(kpc) 


(10^ Mo kpc"^) 


(kpc) 


(10^ Mo) 


(10^ Mo) 


Leol^ 


0.28 ± 0.01 


10.4 ± 1.0 


1.9 


5.3 ± 1.3 


0.27 ± 0.02 


48 


3.1 ± 0.6 


LeoII^ 


0.19 ± 0.01 


7.5 ± 0.6 


0.9 


6.1 ± 1.8 


0.18 ± 0.02 


5.8 


1.1 ± 0.2 


Carina^ 


0.31 ± 0.01 


7.5 ± 0.4 


1.2 


2.1 ± 0.3 


0.32 ± 0.02 


4.3 


1.8 ± 0.2 


Sextans^ 


0.64 ± 0.04 


6.3 ± 1.0 


1.9 


3.5 ± 1.5 


0.65 ± 0.06 


5.0 


2.6 ± 0.8 


Draco"^ 


0.247 ± 0.002 


10.5 ± 0.8 


1.5 


6.9 ± 1.2 


0.24 ± 0.01 


2.6 


2.8 ± 0.4 


Ursa Minor^ 


0.321 ± 0.014 


12.8 ± 1.2 


1.1 


6.6 ± 1.6 


0.28 ± 0.02 


2.9 


5.2 ± 0.9 



Table 7.1: Parameters obtained from mass modeling of six Milky Way dSphs. 
Columns: (1) name of dSph and reference for surface brightness profile used. 
1: [260]; 2: [26T]; 3: [254]; (2) scale-length of Plummer fit to light 
distribution; (3), (4) central velocity dispersion (Tq and 3a lower limit on the 
scale-length Rs of Plummer function fit (equation 17. 3p to velocity dispersion 
profile; (5), (6) central density po and scale-length ro of median Burkert fit 
to density profile from Jeans equations; (7) total stellar mass Mb; (8) total 
mass inside i?83/2 {Rss is the three dimensional radius enclosing 83% of 
the light). Quoted errors indicate la ranges of parameters obtained from 
1000 random realizations of the observed data, but do not include modeling 
uncertainties. The large range of Rs values in each dSph indicates that both 
flat {Rs > lO^kpc) and falling {Rs ~ Ikpc) dispersion profiles are compatible 
with the observations, although in all cases except Ursa Minor, the median 
profile is flat. See j246| for sources of velocity data. 



mer [262] distribution 

where Sq is the central surface density and Rh is the scale-length. Fig. 17.21 
shows the observed surface brightness profiles obtained from the literature 
(see column 9 of table 17.11 for references) and the best-fitting Plummer dis- 
tributions for each dSph in our sample. The median values for i?b obtained 
from the random realizations are given in table 17.11 

We fit each line of sight dispersion profile with a function of the form 

= - — P-^) 
i + (f 

where ctq is the central velocity dispersion and Rs is the scale length of 
the dispersion profile. Table 17.11 also gives all the relevant parameters for 
these fits. Notice that, as in the analysis presented in |246| . the fits to 
the dispersion profiles are merely functional fits to smooth the data and 
reproduce the general shape of the observed dispersion profiles. We observe 
generally flat velocity dispersion profiles, the scale radii Rg being much larger 
than Rh (see column 4 of table [7TT]) . The observation that our measured 
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Figure 7.2: Surface brightness profiles of six Milky Way dSphs. Observed 
profiles are shown as points with error bars. Best-fit Plummer profiles are 
shown as solid curves. In all cases, the Plummer model is a good match to 
the light distribution. 



dispersion profiles for these six dSphs are consistent with being flat, and that 
deviations from this are not statistically significant, is in agreement with the 
more recent (and more extensive) data of [248J for five of these systems (Ursa 
Minor was not included in their sample) in which the dispersion profiles are 
found to remain flat to very large projected radii. Our conclusions in this 
paper would thus be the same if we had used the |248| data. Fig. 17.31 shows 
the observed velocity dispersion profiles, the best-fitting Plummer functions 
to the raw dispersion data and the median Plummer function for each dSph. 

We use the Jeans equations to determine the three-dimensional mass 
profile corresponding to each realization of the light distribution and velocity 
dispersion profile, under the assumptions of spherical symmetry and velocity 
isotropy and fit a Burkert profile to the three dimensional density profiles 
thus obtained. The median values and la ranges of the Burkert parameters 
are presented in table 17.11 Fig. 17.41 presents the density profiles obtained 
from the best-fit surface-brightness and velocity dispersion profile as well as 
the corresponding best-fit Burkert profile. We note that the Burkert profile 
obtained from the median values of po and ro is very similar to the best- 
fit profile shown. Moreover, as a sanity check, in Fig. 17.31 we overplot the 
observed dispersion profile for each dSph with the profile obtained from the 
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Figure 7.3: Line-of-sight velocity dispersion profiles of six Milky Way dSphs. 
Observed profiles are shown as data points with error bars. The solid curves 
show the best-fit Plummer functions (equation 17. 3p to the observed disper- 
sion profiles. The dashed curves show the median Plummer function based 
on the Monte Carlo realizations of the observed data. The dot-dashed line 
shows the dispersion profile obtained using the best-fit Plummer profile to 
the light distribution and the best-fit Burkert model to the halo mass distri- 
bution. 



the best-fit Burkert halo and the best-fitting Plummer light distribution. The 
figure shows that the observed stellar data in each dSph can be reproduced 
by a Plummer distribution of stars embedded in a Burkert halo. 

The aim of the present paper is to investigate whether the scaling laws 
found by |62] for the luminous and DM mass distributions in Spirals are 
compatible with the available dSph kinematics and photometry. To facilitate 
some of this comparison, we must define a stellar length scale for the dSphs 
which plays the same role as the disk scale length i?D in spirals. One way to 
do this is to identify the location of the peak of the Plummer dSph stellar 
spheroid "rotation curve", occurring at 1.4i2b; with the peak, at 2.2. Rd, of 
the stellar Freeman disk RC. Thus, in the dSphs, we associate the Spiral 
length scale Rd with the radius 0.64i?b- However, we note that most of our 
conclusions in this paper do not make use of this length scale. 
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Figure 7.4: Three-dimensional mass density profiles (solid curves) for our 
dSphs obtained using the Jeans equation, assuming spherical symmetry and 
velocity isotropy. The dashed curves show the best-fit Burkert profiles. 



dSph core radii and halo central densities 

Under the assumptions of spherical symmetry and velocity isotropy, the ob- 
servation of a flat velocity dispersion profile implies that the DM mass profile 
inferred via the Jeans equations (e.g. |263| ) is dictated by the distribution 
of the luminous matter: M(r) oc — r d log /9;,(r)/(i log r. If the stellar den- 
sity distribution is cored, e.g. it is represented by a Plummer distribution, 
the DM distribution obtained is also cored with the two core radii being 
proportional. It turns out that the scale length of a Burkert halo fitted to 
the profile obtained in this case is equal to the Plummer radius of the light 
distribution, as can be seen by comparing columns 2 and 6 of table 17.11 

The accurate determination of the size of the DM core radii in dSphs re- 
quires the construction of dynamical models which include velocity anisotropy 
and which can be compared to the full velocity distribution rather than just 
the velocity dispersion as in the Jeans equations. This analysis is beyond 
the scope of the current paper, and will be presented elsewhere in connection 
with a larger velocity data set (Wilkinson et al., in prep.). In the absence 
of such constraints, it is important to consider whether the tq value we use 
are physically meaningful. It is possible that the stellar core radii of dSphs 
may have evolved from their original values since their formation due to 
various processes, both internal (e.g. supernovae) and external (e.g. tidal 
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disturbance). Since our analysis is based on a scale which is essentially the 
present-day core radius of the stellar distribution, we must be cautious in 
drawing conclusions from this about the underlying DM distribution. In 
what follows, therefore, we restrict ourselves to the possible compatibility 
between the values of the the halo parameters as extrapolated at low lumi- 
nosity from spiral galaxy scaling laws and the observed kinematics of dSphs. 
We defer a more robust demonstration of the physical nature of these pa- 
rameters in dSphs (and in particular of their relation to actual DM halo 
parameters) to future work. 

In contrast to the halo core radius ro which is determined from the kine- 
matics and the light distribution, via an assumption on the anisotropy of 
stellar motions, the normalization of the halo density po is constrained by the 
amplitude of the velocity dispersion profile. In our models, po corresponds 
to the mean mass density inside one core radius of the light distribution, and 
is therefore likely to be accurate to better than a factor of three, allowing 
for uncertainty in the velocity anisotropy. 

The assumptions of velocity isotropy and spherical symmetry that have 
allowed us to solve the Jeans equation are supported by additional arguments 
in a two particular dSphs (Ursa Minor and Fornax; see |246| ). Although 
models with larger cores (and appropriate velocity anisotropy profiles) might 
also reproduce the observations, our goal in this paper is to explore whether 
the gross properties of the DM halos around dSphs are consistent with the 
better-determined relations that characterize the ~ 10^ times more massive 
halos around spirals. In this context, it is thus interesting to investigate first 
the consistency of the simplest models. 

7.3 Dark matter properties 

In spiral galaxies, [62] have shown that the DM distribution is closely related 
to that of the luminous matter. Their structural parameters are all corre- 
lated: the mass and the length-scale of the luminous matter correlate with 
similar quantities of the DM ([62])- We start to frame the DM properties in 
galaxies of different luminosity and Hubble Types by analyzing for Spirals 
and dSphs the po vs Rd relationship in Fig. 17.51 which is not affected by an 
anisotropy assumption in a way relevant relevant for our scopes. The data 
are taken from [75] and are in good agreement with those in [TJ [H]. The 
"central" densities of DM halos regularly increase as the size of the stellar 
component decreases. In detail we obtain the intriguing result that although 
dSph halos are much denser, they are found to lie on the extrapolation of 
the spiral relationship. Although the observational evidence for this relation 
is relatively strong, we stress that its physical interpretation is presently 
unknown. 

We continue our comparison between spirals and dSphs by testing the 
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Figure 7.5: Halo central density po versus stellar length scale Kb for spirals 
(solid curve) and dSphs (points) 



dSphs for consistency with the internal halo relationship. We next consider 
how the parameters po and tq for the dSphs, which we discussed in the 
previous section, compare to those of spiral galaxies. In Fig. 17.61 we plot 
Po versus tq, recalling that ro in dSph halos is an assumption-dependent 
quantity and its errorbar does not include the significant uncertainty which 
arises from our lack of knowledge about the velocity anisotropy. Interestingly, 
the figure shows that the extrapolation to higher central densities of the 
Po — ro relation for spirals would predict halos halos for the dSphs which, as 
we have seen, are consistent with the observed kinematics. 

Although the observed data we are using for the dSphs neither require 
cored halos, nor constrain their values in a model-independent way, the ease 
with which a family of dSph halos can be obtained by simple re-scaling of 
larger spiral galaxy halos is intriguing. If confirmed by future data, the exis- 
tence of such a scaling law, spanning three orders of magnitude in each halo 
parameter, would indicate that the physical processes of galaxy formation 
tend to produce DM cores of sizes roughly equal to the stellar cores, in all 
galaxies. This would potentially require a significant revision of our picture 
of galaxy formation: it is difficult to explain the origin of such a scaling 
law in that it relates quantities which do not exist in the standard galaxy 
formation theory (i.e. a core radius and finite central density in the DM 
distribution). 
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7.4 The baryonic-dark matter interplay 

In this section, we investigate the couphng between the distributions of dark 
and luminous matter at the level of the global mass properties. This is easier 
to investigate than the relationship among the structural mass parameters as 
done in the previous section: in both Hubble Types, the global properties are 
less dependent on the modeling assumptions and are less strongly affected 
by observational uncertainties. 

[62] found that the dark and stellar mass inside a reference radius are 
very closely related in Spirals. We now investigate whether this general 
behaviour is also seen in dSphs. We plot the ratio of the stellar mass to 
the halo mass at a radius of Rsa/^ corresponding to the region inside which, 
in spirals and ellipticals, the baryonic matter is always a major component 
of the dynamical mass budget. Moreover, in both spirals and dSphs, the 
baryonic matter inside this radius roughly coincides with the total stellar 
content, the HI content being negligible inside this radius in dwarf spirals 
(see e.g. Figs. 4.13 in [78]). Finally, a convenient coincidence is that this 
radius is approximately the farthest one for which we have kinematic data 
for all objects (dSphs and spirals). 

In Fig. l7.7l we show the well-established result that, in contrast to galaxies 
of other Hubble type (and of much larger stellar mass), dSphs are always 
dominated by DM even in their inner regions. The fraction Mg/Mh sets an 
important physical quantity, namely the percentage of baryonic mass residing 
inside the luminous part of a galaxy and (for dSph also the the percentage 
of baryonic mass tout court, given the general absence of an "external" HI 
component). Bearing in mind that all galaxies are thought to have formed 
with the same initial baryon fraction of roughly 17 per cent [I], the data in 
Fig. 17.71 imply that star formation was very inefficient in processing gas into 
stars in galaxies with stellar masses smaller than W^^Mq. In particular, in 
dSphs we find values for the baryon fraction smaller than 10~^ already at the 
optical radius, that imply even smaller global values (i.e. at their virial radii). 
Fig. 17.71 supports the view that dwarf systems, i.e. objects less massive than 
lO^'^M© (irrespective of their Hubble type) must have experienced massive 
supernova feedback that has strongly limited their star formation efficiency 
(see [77] for a discussion). 

An individuality of the dSphs as compared to other Hubble types is that 
the baryonic fraction at any radius exhibits considerable object to object 
variations, of magnitude about 1 dex. This is several times larger than 
those seen in ellipticals and disk systems [62l [67l 0.291 1264] . and cannot be 
explained by merely the uncertainties in the determination of this quantity 
from the observations. A number of authors have noted that the increase 
of mass to light ratio with decreasing total luminosity seen in the dSphs 
is consistent with a common halo mass scale (interior to ~ 0.6 kpc) but a 
systematically varying baryon fraction |246l 12491 1250] . It is possible that 
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Figure 7.6: Structural halo parameters derived in i) spirals, by means of 
the URC (solid line), and the weak lensing shear (squares), ii) NGC 3741 
(triangle) the darkest spiral in the Local Universe by means of its kinematics, 
iii) Milky Way dSph satellites, by means of their internal stellar kinematics. 
The [?] relation is shown as a dashed line. All these data can be reproduced 
by logpo ~ alogr^ + cost with 0.9 < a < 1.1 



this may arise from environmental effects, perhaps related to their varied 
orbits about the Milky Way, and in particular to their minimum perigalactic 
distances [2651 ISMl Elll EMj • 

7.5 Conclusions 

Dwarf spheroidal galaxies are the lowest luminosity stellar systems which 
show evidence of dynamically significant DM. Moreover, (i) their typical 
stellar masses lie in the range 3x lO^M© to 2x 1O''M0, although the luminous 
masses of some recently discovered objects are as low as lO^M© [269]; (ii) 
the central densities of their DM halos reach almost 10^ times the critical 
density of the universe; (iii) their stellar length scales are of order 0.3 kpc; 
(iv) the DM in these systems typically outweighs the baryonic matter by a 
large factor (from a few tens, up to several hundred). 

Let us stress that all the above quantities are about two orders of mag- 
nitude different from those observed for spiral and elliptical galaxies. There- 
fore, due to these extreme structural properties, an understanding of the 
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Figure 7.7: Ratio of the stellar mass to the halo mass at versus the 

total mass inside Rssf^ for the dSphs (points). The solid curve represents 
the relation obtained for spires. 

formation of dSphs is crucial for the development of a complete picture of 
galaxy formation. 

The main result of this paper is the finding that these galaxies, despite 
being very separate in their physical properties from spirals and ellipticals 
and having a large individual scatter in their baryonic properties, exhibit 
kinematics which could be consistent with the presence of DM halos which 
are essentially scaled-down versions of those found in galaxies of much higher 
mass and different Hubble type. We have shown that a Burkert halo density 
profile can reproduce the available kinematic data for the dSphs. We find 
that the derived central densities and the stellar core radii are consistent 
with the extrapolation of the relationship between these quantities seen in 
spiral galaxies. In addition, we have shown that if we extrapolate the rela- 
tion between halo central density and DM core radius previously found in 
ellipticals and spirals, the halo parameters expected for the dSphs would be 
consistent with their observed kinematics. 

This potential consistency is intriguing, and could point to a common 
physical process responsible for the formation of cores in galactic halos of all 
sizes, or to a strong coupling between the DM and luminous matter in dSphs. 
If confirmed, this would suggest a Grand Picture for galaxy formation in 
which in galaxies of all Hubble Types, the DM is "aware" of the length scale 
of the luminous matter and vice versa. It is worth noting that a potential 
connection between spiral galaxies and dSphs does not appear as natural as 
one between dSphs and other hot, spheroidal systems. For example, while 
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the sizes of spiral galaxies are presumably fixed by the angular momentum of 
the gas from which they form, most of the present-day dSphs show no signs 
of rotation ( |259| have recently found evidence of rotation in the Sculptor 
dSph). However, [270j have proposed a formation scenario for dSphs in which 
they are initially low-mass disk galaxies that are subsequently transformed 
into spheroids by tidal interaction with the Milky Way. More recently, such 
models have been shown to provide reasonable models for the properties 
of the Fornax |271| and Leol |272| dSphs. If the halos of dSphs do indeed 
follow the scaling laws defined by more massive disk galaxies, this could lend 
indirect support to evolutionary histories of this kind. 

We also find evidence that the depletion of primordial gas through su- 
pernova feedback has proceeded in a similar manner across all Hubble types, 
with the resulting luminous to DM ratio depending mostly on the depth 
of the gravitational potential. In the dSph potentials, which correspond to 
a virial temperature of order lO^K, we find a depletion by a factor of one 
hundred at Rss, and a factor which may reach and exceed 10'^ at the virial 
radius. 

As we have emphasized throughout this work, further dynamical analysis 
is needed in the dSphs to show directly that they possess DM cores and, if 
so, to constrain their core radii. Nevertheless, it interesting to speculate on 
the possible implications of these scaling laws for our understanding of DM. 
Warm DM has been invoked as a potential solution to the over-prediction of 
substructure by ACDM simulations, and to the cusp-core issue (e.g. |164| ). 
However, the existence of scaling relations between the central density and 
core radius over three orders of magnitude in both quantities would rule 
out this explanation, unless the warm DM spectrum is extremely fine-tuned. 
Further, such DM relations cannot arise due to either self-annihilation or de- 
cay of DM which would predict a narrow range in po and no clear correlation 
of the latter with the core radius. 

|166| argued that the phase-space densities of DM halos suggested that 
warm DM (either collisional or collisionless) could not be the cause of cores 
in galaxy halos on all scales. These authors suggested a dynamical origin 
for the cores of larger galaxies. A universal scaling relation suggesting that 
any core formation process has to proceed with approximately comparable 
efficiency across three orders of magnitude in scale, would render dynamical 
core formation scenarios (e.g. angular momentum transfer from the baryons 
to the halo, expulsion of baryonic matter by supernovae, or spiralling binary 
black holes at the centre of the galaxy, etc.) more difficult to envisage. We 
can speculate that a physical property of DM which has the potential to ex- 
plain the origin of the observed trends among the structural DM parameters 
in primordial NFW halos would be a self-interaction with an appropriate 
velocity- dependent cross-section. Alternatively, some currently unknown in- 
teraction between DM particles and baryonic matter or photons may be 
required to explain core formation at the galactic scale. 
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Clearly, direct kinematic evidence for or against the presence of cores in 
dSph halos is now required to resolve the situation. If cores are detected and 
are found to have parameters consistent with those discussed in this paper, 
this will provide important information about the properties of the DM of 
which they are composed. On the other hand, if it turns out that cored 
halos are not a general feature of dSphs (restricted perhaps to the cases of 
the Ursa Minor and Fornax dSphs which require cores to allow survival of 
their internal substructure), the similarity of the apparent interplay between 
dark luminous matter in dSphs and spirals, as suggested by Fig. 17.51 would 
remain an intriguing observation. 



Chapter 8 



A constant Dark Matter Halo 
Surface density in Galaxies 

In the same line of the previous chapter I investigate further the DM halo 
properties. In particular I discuss our work where we confirm and extend an 
earlier claim by Spano et al. 2008 [H] that the central surface density /ioD of 
galaxy DM halos is nearly constant, independent of galaxy luminosity. Based 
on the co-added RCs of ~ 1000 spiral galaxies, mass models of individual 
dwarf irregular and spiral galaxies with high-quality RCs, and the galaxy- 
galaxy weak lensing signals from a sample of spiral and elliptical galaxies, we 
find that log^oD = 2.05 ± 0.15, in units of M© pc~^. We also show that the 
observed kinematics of Local Group dwarf spheroidal galaxies are consistent 
with this value. Our results are obtained for galactic systems spanning a 
wide range in magnitude, belonging to different Hubble Types, and whose 
mass profiles have been determined by independent modeling methods. The 
constancy of /Uqd is in sharp contrast to the variation, by several orders 
of magnitude, of the halo density and stellar surface density in the same 
objects. 

8.1 Introduction 

It has been known for several decades that the kinematics of disk galaxies 
exhibit a mass discrepancy: in their outermost optical regions the circular 
velocity profile cannot be explained by the ordinary stellar or gaseous matter. 
This is usually solved by adding an extra mass component, the DM halo. 
RCs have been used to assess the existence, the amount and the distribution 
of this dark component (e.g. [HIE^)- Recent debate in the literature has 
focused on the "cuspiness" of the DM density profile in the centres of galaxy 
halos that emerges in CDM simulations of structure formation [3j 11461 11641 
[273] but is not seen in observed data (e.g. [6l[l0l[IIIl[Il2l[l33l[l79]), as well 
as on the various systematics of the DM distribution (see [75]). A significant 
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contribution to this debate was recently made by Spano et al. 2008 [H], who 
fitted the RCs of 36 spiral galaxies using a mass model involving a cored dark 
sphere of density 



where pQ is the central density and tq is the core radius. The authors found 
that the quantity ^uqd = PqI'o, proportional to the central halo surface density 
Ti{R) = 2 p{R, z)dz, is independent of the galaxy blue magnitude: 



log(yUoD/M0pc~^) = 2.2 ± 0.25 or pqd = 150^71^° Mqpc"^. (8.2) 



For the sake of completeness, we note that a constant pqd of about 100 Mqpc" 
but with a much larger r.m.s (0.4 dex), was found in the earlier work of |274| . 
for a sample of 50 spiral and dwarf galaxies. 

In this work, we will investigate the constancy of poD found in [H] for ob- 
jects whose central densities and core radii vary by 1-2 orders of magnitude. 
We aim to confirm or rule out this property by investigating independent 
samples of galaxies that include a large number of objects of different Hub- 
ble Type and magnitude and whose halo properties have been estimated 
using different and independent methods of mass modeling. Given the wide- 
ranging nature of the data and models we include, a positive result arising 
in this study would be difficult to dismiss as a coincidence. 

In this work, we make use of data from galactic systems spanning wide 
ranges in luminosity and Hubble Type. Moreover, their mass distributions 
are modeled by means of different techniques. In particular, our results 
are obtained from mass models of: (a) a large sample of Spiral galaxies, 
analyzed by means of their URC; (b) the darkest Spiral in the local Universe, 
studied through its kinematics; (c) a large sample of Spiral and Elliptical 
galaxies, for which weak-lensing shear measurements are available. We also 
compare the value of pqd obtained from these luminous galaxies with the halo 
parameters consistent with the kinematics of six dwarf spheroidal satellite 
galaxies of the Milky Way for which extensive stellar kinematic data sets 
are available. We note that with the exception of the weak lensing results 
which are presented in this work, the values of po and ro (and their relative 
uncertainties) that we use to compute pqd and then to investigate Eq. 18.21 
are obtained and discussed in previous works. We will refer interested readers 
to those publications for details of the data and models. 

In Section 2, we compute the quantity pqd for different families of galax- 
ies and compare it with the [IT] result. A discussion of our result is given in 
Section 3. 



p{r) 




3/2' 



(8.1) 



8.2. THE pqRq VS MAGNITUDE RELATIONSHIP 



105 



8.2 The pqTo vs magnitude relationship 

In this work, we assume that the DM halo in each galaxy follows the Burkert 
profile ^08j: 

= , ^ u'L 2v (8-3) 
(r + ro) [r^ + rg) 

This profile, when combined with the appropriate baryonic gaseous and stel- 
lar components, has been found to generally reproduce very well the available 
kinematics of disk systems out to 6 Rd (IHEIE]; see [TU] for the case of the 
most extended RC). 

The possible existence of a constant central surface density of DM for 
all galaxies does not depend on which specific (cored) density distribution 
we assume for the DM, whether we adopt any of the following: Spano et 
al. (2008; labeled as S hereafter), Donato et al. (2004; D) [7] or the 
present one (B). Since different cored mass models provide equally good 
fits to the same kinematical data sets (e.g. [6]), with all of them (presum- 
ably) describing the true, underlying halo mass profile M^, the relations 
Mh{r,B) = Mh{r,S) = Mh{r,D) = Mh{r,true) must hold, to within ob- 
servational uncertainties. This enables us to derive proportionality factors 
between the corresponding parameters of the different cored profiles. These 
can easily be computed: log fiQD{D) = log fiQ£){B) + 0.1 = log^oD('S') + 0.3 
showing that the correction terms needed to compare different profiles are 
quite negligible for each specific profile, relative to the observed object-to- 
object variance of /xqd at a fixed magnitude. 

Let us consider the case in which the halos around galaxies are cusped 
rather than cored, as predicted by cosmological simulations of structure for- 
mation (e.g. [3l I275| ). then does the use of the Burkert profile introduce a 
bias into the results we obtain? We first remind that this that possibility is 
unlikely in view of the many cases in which the NFW profile fails to fit the 
observed spiral kinematics. However, in any case, in the range 0.2R]:) — R^r, 
the Burkert profile (with a small value for the core radius and a appropriate 
value for the "central density"), can mimic quite well the velocity profile of 
a NFW halo with a standard value of the concentration parameter. The 
Burkert profile is therefore an empirical one able to "measure" the level of 
cuspiness of the underlying DM density distribution. As general result, with 
the same number of free parameters (i.e. a length scale and a density scale) 
the Burkert profile is able to fit all available kinematical data within the 
observational uncertainties; moreover, differently from NFW mass model- 
ing the present one is able to estimate very properly the disk mass, which 
turns out in agreement with the expectations from stellar population syn- 
thesis models (e.g. [6l[IU[80], see also |234| ). Thus, the halo parameters we 
use in this work are suitable and unbiased measures of the spirals physical 
properties. 
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Immediate, though indirect, support for the [11] claim comes from the 
results of [7]. In Fig. 1 we plot ^qd as a function of the stellar exponential 
scale-length Rd for the sample of 25 disk systems (Spirals and LSB) analyzed 
by We see that the derived values for ^oz) are almost constant, although 
Rd varies by more than one order of magnitude. In addition, there is no 
obvious difference between the results from High Surface Brightness (HSB) 
galaxies and Low Surface Brightness (LSB) galaxies. This result is good 
agreement with Eq. 18.21 However, it is important to note that the two 
samples are similar, with five objects in common, and the analysis employed 
is essentially the same. 
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Figure 8.1: The central halo surface density poTq as a function of disk scale- 
length Rd for the Donato et al. (2004) sample of galaxies. Open and filled 
circles refer to LSB and HSB galaxies, respectively. The solid line is our best 
fit to the data. 

We now calculate the central surface density /xqd for the family of Spirals 
by means of their URC. This curve, on average, reproduces well [MlCS] the 
RCs of individual objects out to their virial radii Ryir (the radius at which 
the halo mass is 100 times the background mass). The URC is built from 
(a) the co-added kinematical data of a large number of Spirals (|62j; see also 
|76| ) and (b) the disk mass versus halo virial mass relationship found by |77] 
and it leads, for objects of given luminosity (or disk mass), to specific values 
of po and tq (see equations 6a, 7 and 10 of [75] for details). The solid line in 
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Figure 8.2: poro in units of Mqpc~^ as a function of galaxy magnitude for 
different galaxies and Hubble Types. The original Spano et. al. (2008) data 
(empty small red circles) are plotted alongside those for Spirals obtained by 
the URC (solid line). The full (green) circle corresponds to the the dwarf 
galaxy N3741, full (black) squares to Spirals and Ellipticals obtained by weak 
lensing and the (pink) triangles to the dSphs obtained by their kinematics. 



Fig. 2 shows the resulting hod as a function of galaxy magnitude. Because 
the URC is derived from co-added RCs, the particularities of individual 
galaxy curves (e.g. observational errors or non-axisymmetric motions due to 
bars or spiral structure) are averaged out. The URC therefore allows us to 
trace the general form of the gravitational potential of Spirals over their full 
luminosity range. A natural concern is that the values of the halo parameters 
we obtain are biased by the smoothing process itself. However, the values 
of Pod obtained from detailed mass modeling of 36 RC of spirals by [TT] 
shown in Fig. 2 as open circles (the pqd's for the 25 mass models in [7] 
not reported here are in very good agreement with the latter) are consistent 
with those obtained from the the URC, all suggesting that these various mass 
modeling it is returning physically meaningful values of physically meaningful 
mass parameters. More in detail, the URC provides, for Spirals of a given 
luminosity, a reliable estimate of their average value of poD, although not 
of their cosmic variance around it. In the estimation of the latter quantity, 
the detailed studies of individual objects such as those of [H] and [7| are 
indispensable to provide us with the needed quantity, that results negligible 
for the present aim. 

In this work, we estimate the values of pqd from the DM structural pa- 
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Figure 8.3: Tangential shear measurements from |276| as a function of pro- 
jected distance from the lens in five B-band luminosity bins. In this sample, 
the lenses are at a mean redshift z~0.32 and the background sources are, in 
practice, at 2: = oo. The solid (dashed) magenta line indicates the Burkert 
(NFW) model fit to the data. 

rameters obtained in a third different way: from analyzing the galaxy-galaxy 
weak-lensing signals for a large sample of Spiral and Elliptical galaxies. The 
details are presented as follows. 

Recent developments in weak gravitational lensing have made it possi- 
ble to probe the ensemble- averaged mass distribution around galaxies out 
to large projected distances. These new data provide crucial information, 
complementary to that obtained from kinematics. The tidal gravitational 
field of the DM halos generates weak-lensing signals, by introducing small 
coherent distortions in the images of distant background galaxies, which can 
be detected in current large imaging surveys. We can measure, from the 
centre of the lenses out to large distances (much greater than the distances 
probed by the kinematic measurements), the azimuthal-averaged tangential 
shear 7t 

S(m-S(i?) , 
< 7t >^ ^ ' ^ \ (8.4) 

where S(i?) = 2 p{R, z)dz is the projected mass density of the object dis- 
torting the galaxy image, at projected radius R and S(i?) = J^^ xT,{x)dx 
is the mean projected mass density interior to the radius R. The critical 
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density Sc is given by 



A-kG a As 



(8.5) 



where A and D\ are the distances from the observer to the source and 
lens, respectively, and As is the source-lens distance. The above relations 
directly relate observed signals with the underlying DM halo density. For our 
analysis we use the weak lensing measurements from |276| available out to a 
projected source-lens distance of 530 kpc. The sample, which contains about 
10^ isolated objects and spans the whole luminosity range of Spirals, is split 
into 5 luminosity bins of magnitudes given in table 18.11 The most luminous 
bin is likely dominated by the biggest Ellipticals. By adopting a density 
profile, we model 7t (see Fig. 18. 3p and obtain the structural free parameters 
po and ro by means of standard best-fitting techniques. The Burkert profile 
given by equation 18.71 provides an excellent fit to the tangential shear (see 
Fig. 18.31 and table 18. ll) . The NFW density profile provides a less satisfactory 
fit to the gravitational shear around the most luminous objects (Fig. 18.31 see 
also Fig. 6 of |276| . Notice that at fainter luminosities {Mb > —20.1) the 
signal-to-noise is too low to discriminate between mass models, so that while 
the Burkert profile remains a working assumption, NFW profiles cannot be 
excluded. Assuming the Burkert halo profile we plot the resulting /Uqd values 
in Fig. 2 as solid squares. 
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Table 8.1: Structural parameters and goodness of fit for a Burkert profile to 
the weak lensing signal of [276j . 



The nearby dwarf galaxy NGC 3741 {Mb = —13.1) is a very interesting 
case: it represents the very numerous dwarf disk objects which are DM dom- 
inated down to one disk length-scale or less and in which the HI gaseous disk 
is the main baryonic component. In addition, this specific galaxy has an ex- 
tremely extended and very symmetric HI disk, which allowed [TO] to carefully 
trace the RC and therefore its gravitational potential out to unprecedented 
distances relative to the extent of the optical disk. The data probe to radii 
of 7 kpc (equivalent to 42 B-band exponential scale lengths), and have sev- 
eral independent points within the estimated halo core radius. The RC was 
decomposed into its stellar, gaseous and dark (Burkert) halo components, 
yielding a very good fit |180j : the corresponding //qd is plotted in Fig. 2 as 
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a filled circle. The relatively large error-bar is due to uncertainties in the 
distance. 

At the level of 0.2 dex, no large differences emerges between the values 
of //QD estimated in different way or referring to a Spiral or an Elliptical 
population. It thus appears that the central surface density of DM halos 
assumes a nearly constant value with respect to galaxy luminosity, over a 
range of at least nine magnitudes. 

The Milky Way satellite dwarf spheroidal (dSph) galaxies are the smallest 
and most DM dominated systems known in the universe (see e.g. [2461 1249j 
and references therein). Their low HI gas content is another property that 
sets them apart as a galaxy class (e.g. |277| ). In a recent study of six dSphs, 
|246| showed that, assuming spherical symmetry and velocity isotropy, the 
stellar kinematics and photometry of dSphs are consistent with their occu- 
pying cored DM halos. [278] subsequently showed that, for the same sim- 
plifying assumptions, cored Burkert profiles are able to reproduce the dSph 
kinematic observations. Our current lack of knowledge about the anisotropy 
of the velocity distribution means that the density profiles of dSphs are not 
uniquely constrained by the data, and both cored and cusped models can re- 
produce the data in most dSphs |246ll250ll259| . Bearing this caveat in mind, 
it is nevertheless interesting to compare the value of /ioD from our spiral and 
elliptical galaxy samples with the values obtained from the halo parameters 
which |278| showed to be consistent with the dSph kinematics. These are 
plotted in Fig. 2 as triangles. Note that the errorbars shown reflect only 
the statistical errors in the estimation of the parameters from the observed 
data, and do not account for any modeling uncertainties. We emphases that 
the relatively small range of both halo density and core radius found for the 
dSphs means that the current data in these galaxies would be consistent 
with the approximate constancy of any product of po and rg. In particular, 
it has been noted that all the dSph data are consistent with their occupying 
halos which contain roughly equal masses interior to about 0.6-1.0 kpc (i.e. 
po^^o ^ constant: |246l 1279] . Nevertheless, the consistency of the dSph data 
with the value of poD suggests that the relation pqVq ^ constant may extend 
to fainter systems, and thus be valid over a range of fourteen magnitudes in 
luminosity. 

8.3 Discussion and Conclusions 

We have compiled data on the DM halo mass distribution in many galactic 
systems of different Hubble Type (including Dwarfs disk galaxies. Spirals, 
Ellipticals) spanning a luminous range of about —8 < Mb < —22 and a 
gaseous-to-stellar mass fraction of many orders of magnitude. The mass 
modeling of such objects has been carried out using different and independent 
methods. The halos are all well reproduced by a cored profile with two 
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structural parameters: a central halo density po and a core radius tq, whose 
respective values range over several orders of magnitude: 6 x lO^^^g/cm^ < 
Po < 10~^^g/cm^ and 0.3 kpc < tq < 30 kpc. In spite of dealing with 
galaxies with such different physical properties, we have found that their 
central DM surface density hqd = PoTq remains almost constant: 



independent of galaxy luminosity. In addition, we have compared this value 
of Pod with the kinematics of dSphs and found that these are also be con- 
sistent. 

Our results support the pioneering analysis by [H], based on a sample of 
36 spiral galaxies, in which they found a nearly constant halo surface den- 
sity around 150 MqPc"^ independent of galaxy luminosity. In addition to 
investigating many more objects across more Hubble- types and a much wider 
luminosity range, we have obtained the halo surface density poD both from 
individual galaxy kinematics and from co-added kinematical/shear measure- 
ments. The approximate constancy of pon is in stark contrast to the stellar 
central surface density in galaxies of different Hubble Type and magnitudes 
which shows large variations (see the relevant works cited above for de- 
tails). In Spirals, it ranges between SOOMqpc"^ at about Mb = —22.5 to 
~ 5OM0pc~^ at Mb = —17, in dSph it probably does not reach IMqpc"^, 
while in Ellipticals it easily exceeds 10000 Mqpc"^, with large variations 
with luminosity and object-to-object. 

It is important to consider how the approximate constancy of poD with 
Mb is related to the correlation between tq and po, 



which has been claimed in Spiral galaxies [108j . First, the former relationship 
(Eq. 18. 2p links two very different physical properties of galaxies (i.e. the 
central DM surface density and galaxy magnitude), while the latter (Eq. 
18.7(1 relates two "internal" DM halo structural parameters. Further, let us 
stress that A could be near, but maybe relevantly not coinciding with -1, see 
[4l [6l iTl [751 1108] and even show some non (log) linearities (see [6]), but still 
the quantity rg po could be found constant, within a factor 2, over several 
orders of magnitudes. The study of Eq. 18.61 and Eq. 18.71 must therefore 
proceed separately. 

The evidence that the DM halo central surface density PoTq, over at least 
nine (and possibly up to fourteen) galaxy magnitudes and across several 
Hubble types, remains constant to within less than a factor of two, suggests 
that Pod may be an important physical quantity in the DM distribution 
of galaxies. This is a surprising finding, as it is difficult to envisage how 
such a relation can be maintained across galaxies which range from DM- 
dominated to baryon-dominated in the inner regions. In addition, these 



POD = 110ti[! M0pc"2 



(8.6) 



log ro = A log po + C 



(8.7) 
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galaxies have experienced significantly different evolutionary histories (e.g. 
numbers of mergers, significance of baryon cooling, stellar feedback, etc.). 
Further investigation is clearly required in order to verify and interpret this 
relation. 



Chapter 9 

Conclusions 



Although CDM model is able to correctly describe observations made on the 
largest cosmological scales down to roughly those of galactic scales, and from 
the early Universe to the present epoch, on subgalactic scales it predicts that 
there should be more DM than is detected gravitationally. 

There are several suggestions that could account for the lack of the cusps, 
both from fundamental physics and through astrophysical processes. These 
suggestions make definite predictions of other observables that could be used 
to test the variant properties of DM. In this Thesis I show how these tests can 
be performed with the analysis of the RCs with a proper sample of spirals. 

As gravity is by far the dominant interaction at cosmological scales and 
the force governing the evolution of the universe, another perspective to the 
current picture of the evolution and the matter content of the Universe arises: 
the description of the gravitational interaction at the relevant scales may be 
not sufficiently adequate and a modification of gravity could answer the cos- 
mological and astrophysical riddles. Conceivable alternatives are numerous 
and each of them produces distinctive modifications on small scales that can 
be tested through improved astronomical observations and numerical sim- 
ulations. In my Thesis I explore how such modifications may well account 
for the phenomenon of the RCs. Of course such a solution pays the price of 
renouncing the great success of the actual theory of structure formation and 
evolution envisaging a pure baryonic scenario. 

In my Thesis I also investigate the extension of the well-known scaling re- 
lations of DM halo properties on a large range in galaxy luminosity, including 
the latest observations of the Milky Way satellites. Within some assumption 
I show that the halos of the faintest objects have properties scaled down with 
respect of the bigger ones. This possible consistency could well point to a 
common physical process on the formation of galactic halo cores of all sizes. 
I also find a surprising evidence of a DM constant halo central surface den- 
sity over a large range in galaxy magnitude, suggesting that it could be an 
important physical quantity in the DM distribution of galaxies, even though 
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these galaxies have different evolutionary histories. Further investigations 
are necessary to better interpret these results. 

I conclude that properties of dark matter are written on the kinematical 
features of the luminous matter and that their detailed study can give the 
right glasses to decipher its nature. 
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